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Abstract 

Several A'^-body problems in ordinary (S-dimensional) space are introduced which are 
characterized by Newtonian equations of motion ("acceleration equal force;" in most 
^ I cases, the forces are velocity-dependent) and are amenable to exact treatment ("solv- 

■ able" and/or "integrable" and/or "linearizable"). These equations of motion are al- 
I ways rotation- invariant, and sometimes translation-invariant as well. In many cases 
. they are Hamiltonian, but the discussion of this aspect is postponed to a subsequent 

<^ \ paper. We consider "few-body problems" (with, say, Af=l,2,3,4,6,8,12,16,...) as well 

. as "many-body problems" (A'^ an arbitrary positive integer). The main focus of this 

paper is on various techniques to uncover such A'^-body problems. We do not discuss 
i the detailed behavior of the solutions of all these problems, but we do identify several 

models whose motions are completely periodic or multiply periodic, and we exhibit in 
(-H . rather explicit form the solutions in some cases. 

> 



X 



1 Introduction 

In this paper we exhibit several solvable and/or integrable and/or linearizable dynam- 
ical systems which are naturally interpretable as NAiody problems in ordinary (three- 
5^ \ dimensional) space, inasmuch as they are characterized by equations of motion of Newto- 

nian type ("acceleration equal force"), 

rn = fn{rm,rjn]m = l,...,N^ . (1.1) 

Here and throughout, superimposed arrows denote 3-dimensional vectors, and super- 
imposed dots denote differentiation with respect to the independent variable t ("time"). 
Of course the 3-vector rvi {t) identifies the position of the n-th particle at time t. Unless 
otherwise indicated, the index n, as well as most other indices, runs from 1 to N . In some 
cases is a small integer, say = 1, 2, 3, 4, 6, 8, 12, 16 ("few-body problems"); in others, 
it is an arbitrary positive integer ("many-body problems"). 
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Our main purpose in this paper (the first of a series) is to introduce techniques whereby 
many such models can be manufactured. Only some of these models are actually exhibited 
herein; and only in some cases do we discuss in any detail the solutions of their equations 
of motion, namely the actual motions entailed by the A'^-body problem under consider- 
ation. Yet we trust our treatment demonstrates the interest of these models. It should 
be emphasized in this respect that, until recently, essentially all known solvable and/or 
integrablc and/or lincarizabic #-body models were embedded in one-dimensional space 
(except for the important but essentially trivial case of a collection of harmonic oscilla- 
tors, characterized by linear equations of motion). Only recently some solvable and/or 
integrable many-body problems in two dimensional space were introduced [1], as well as 
some solvable few-body problems in 2- and 3-dimensional space [2], by techniques com- 
pletely different from those introduced in this paper. The dynamical behaviors exhibited 
by these models are much richer than those characterizing models in one-dimensional 
space, where scattering processes are, to say the least, not very spectacular, and particles 
cannot get past each other without colliding (which is instead not the case already in two- 
dimensional space). Indeed, a number of interesting phenomena are exhibited by solvable 
models in the plane: nontrivial scattering, nontrivial periodic motions, limit cycles, etc. 
[1] . Three-dimensional space allows for an even richer gamut of motions; moreover, we 
happen to live in three-dimensional space, and wc are therefore specially interested in 
models that evolve in the same environment. Indeed, wc suggest that some of the models 
treated herein are suitable to test the reliability of numerical codes for the integration of 
many-body motions in three-dimensional space, such as those used in "realistic" molecular 
dynamics calculations. Some of these models are also likely to have applicative relevance, 
but we do not elaborate on this aspect in the present paper. 

These remarks underscore the relevance of the results reported herein, and justify our 
exclusive interest, in this paper, on the 3-dimensional case (although it will be clear in the 
following that certain results are easily extendable, for instance, to 2- and 4-dimensional 
spaces). But let us also emphasize that the results reported herein consist mainly of 
appropriate reinterpretations of more or less known findings concerning treatable models 
for the time-evolution of matrices as well as appropriate parameterizations of matrices in 
terms of 3-vectors, the combination of these two techniques to yield A^-body problems in 
3-dimensional space being, as explained below, the main new idea underlying the approach 
introduced in this paper. 

We deem that evolution equations of type (1.1) qualify as Newtonian equations of 
motion for an iV-body problem in 3-dimensional space iff they are rotation-invariant; it is 
otherwise trivially easy to manufacture solvable and/or integrable A^-body problems in d- 
dimensional space, with arbitrary d, by an appropriate reinterpretation of known solvable 
and/or integrable M-body problems in one- dimensional space, with M = Nd. 

Indeed, the equations of motion of type (1.1) we consider below are all rotation- 
invariant, namely the "forces" /„ can generally be written as follows: 

N N 

fn ~ ^ ^ [fm^nm "I" fm'pnm] ~l" ^ ^ i^mi ^ fm2fnm\m2 
m=l mi,m2=l 

"I" fmi ^ fm2 'Pnmim2 "I" fmi ^ fm2 'Pnmim2] j 

with the scalar functions (p, (p and (p depending on the A'' {2N -\- 1) (scalar) products fj ■ 



(1.2) 
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f^j ■ f^k, ■ and on the (A^ — 1) (3A^ — 2) /2 (pseudoscalar) triple products fj ■ 
rk /\rm, Vj ■ rk ^fm-, rj ■ fk ^^m, rj ■ fk /\rm- This is the most general rotation-invariant 
force one can write in 3-dimensional space. In fact generally (see below) only a few of 
the scalar functions Lp,ip or cp are actually present (most vanish identically), and moreover 
those few which arc present only depend on a few of the scalar and pseudoscalar products 
listed above (and generally these scalar functions only depend on few of their subscripts, 
see below). 

Some of the models given below are moreover translation-invariant, namely they are 
characterized by forces fn which remain invariant under the translation 

rn^rn + ro , (1.3) 

with rb an arbitrary constant (time-independent) 3-vector. 

Some of the models treated below involve scalar quantities evolving in time, in addition 
to three- vectors. 

Several of the models considered below are Hamiltonian, but we defer (with a few 
exceptions) the discussion of this aspect to a subsequent paper of this series. 

Most of the models considered below feature "velocity-dependent" forces, namely the 
forces /„ in the right-hand-side of (1.1) do depend on (at least some of) the velocities r*„j ; 
this dependence is generally linear or quadratic. The forces depend moreover on (at least 
some of) the "particle coordinates" f^n ; this dependence is generally nonlinear, in most 
cases strongly so (rational). 

Some of the models considered below feature forces that arc homogeneous functions of 
degree one of their arguments, so that (1.1) is invariant under the "scale transformation" 

CVn , (1-4) 

with c an arbitrary (of course scalar) constant. 

Several of the models considered below feature only few-body forces. Particularly in- 
teresting are the cases with two-body forces only, 

TV 

fn ~ ^ ^ fnm (^^njfn'Tfmjfm^ ■ (-'-•^) 
m=l 

We never consider "trivial" cases with "triangular" forces which vanish for m < n oi 
for m > n . Indeed in some cases the two-body forces fnm in the right-hand-side of (1.5) 
are independent of the indices n and m, except for the property of acting only among 
"nearest neighbours" (which of course need not be physically close to each other as they 
move), say 

fnm (^fm fn'i fmj 1"m^ ~ ^m,nf i^n] ^n^ "I" ^m,n+lf^'^^ (^^m fn'i '"n-l-lj ^n+l^ 

+ Sm,n-lf~^ (r'n,rn;rn-l,fn-l^ ■ (1.6) 

(In these cases with "nearest-neighbour" forces one should also add "boundary conditions" 
specifying the behavior of the "first" and "last" particles; we will generally omit to do so 
in this paper, since we stop short of discussing in detail the behavior of the various models 
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we introduce). And we also have cases with two-body forces altogether independent of the 
particle indices, 

including cases where / has a very simple form (see below). 

In any case, for the convenience of the hasty browser we provide, in the following Sec- 
tion ^, a representative list of the solvable and/or integrable and/or linearizable A/^-body 
problems treated in this paper. We consider convenient to make the heuristic (useful if not 
quite precise) distinction among "solvable," "integrable" and "linearizable" models, as fol- 
lows. Solvable models are characterized by the availability of a technique of solution which 
requires purely algebraic operations (such as inverting and diagonalizing finite matrices), 
possibly solving known linear, possibly nonautonomous, ODEs in terms of known special 
functions (say, of hypergeometric type), and perhaps the inversion of known functions (as 
in the standard solution by quadratures). Integrable models are those for which a "Lax 
pair" approach is available, yielding an adequate supply of constants of the motion. As 
a rule these latter models are also solvable, but generally this requires more labor; the 
solvable models are generally also integrable. In the Hamiltonian cases, all these models 
are generally integrable in the Liouville sense. We moreover introduce the notion of lin- 
earizable problems whose solution entails, in addition to the algebraic operations required 
to treat solvable or integrable models, the solution of linear, generally nonautonomous, 
ODEs, which, in spite of their being generally rather simple, might indeed give rise to 
quite complicated (chaotic?) behaviors. In the Hamiltonian cases these latter many-body 
problems need not be integrable in the Liouville sense, although the linearity of the equa- 
tions to be finally solved entails the possibility to introduce constants of the motion via the 
superposition principle, which guarantees that the general solution can be represented as 
a linear combination with constant coefficients of an appropriate set of specific solutions. 
In any case, a linearizable many-body problem is certainly much easier to treat than the 
generic many-body problem, inasmuch as its solution can be reduced to solving a linear 
first-order matrix ODE (indeed, in many cases, a single linear second-order scalar ODE 
- albeit a nonautonomous one - see below). 

Clearly these three categories of problems - solvable, integrable, linearizable - are or- 
dered in terms of increasing difficulty (to solve them), so that, as indeed our language 
(including the title of this paper) indicates, problems belonging to a lower category gen- 
erally also belong to the following ones. 

As much as we tried above to delineate clearly the distinction between solvable, inte- 
grable and linearizable models, we are well aware that this categorization is imprecise; for 
one thing, we left vague what "solving" a problem means: finding the general solution? 
solving the initial-value problem? for which class of initial data? and what about bound- 
ary conditions (which in some cases, see above and below, are essential to completely 
define the problem)? Moreover, as indicated above, we are well aware that the boundaries 
among the three categories we have introduced lack complete cogency. Yet this distinction 
provides an useful means to convey synthetically significant information on the status of 
the different models we introduce, a status whose more precise definition will then require 
a detailed case-by-case analysis of the problem under consideration, which clearly exceeds 
the scope of this presentation (lest its size approximate a book rather than a paper). 



(1.7) 
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Our starting point to obtain solvable and/or integrable and/or linearizable A'^-body 
problems in ordinary (3- dimensional) space are solvable and/or integrable and/or lin- 
earizable evolution equations for matrices. A survey of such models is provided in Section 
^ (which might appeal to some readers in its own right), and a number of matrix represen- 
tations in terms of 3-vectors suitable to generate rotation-invariant A^-body problems in 
3-dimensional space are reviewed in Section |^. The material of these two sections is then 
used in Section ^ to derive several solvable and/or integrable and/or linearizable A^-body 
problems in 3-dimensional space, including, in some cases, their explicit solutions and/or 
some information about their behaviors (for instance, the identification of cases with com- 
pletely periodic solutions). As mentioned above, a representative list of these problems is 
also reported in the following Section |2|: we deem the advantage thereby gained in terms 
of clarity of presentation, especially for the hasty browser, to overcompensate for the in- 
evitable repetitiveness this entails. A final Section ^ outlines future researches suggested 
by the findings reported in this paper. To avoid the flow of the discourse being too clut- 
tered by detailed computations some of these are relegated to a number of appendices, as 
well as some flndings that do not quite fit in the mainstream of our presentation. 

2 List of solvable and / or integrable and / or linearizable few- 
body and many-body problems 

In this section we display, for the convenience of the hasty browser, hence with minimal 
commentary, a representative (namely, far from complete) list of the A^-body problems 
treated in this paper. To maximize the user-friendly character of this presentation, we 
separate few-body from many-body systems, and we indicate in each case whether the 
equations of motion presented are solvable, integrable or linearizable (for the explanation 
of this distinction see the end of the introductory Section |^). Some of the many-body 



problems exhibited in Section 2^ are direct generalizations of corresponding few-body 
problems listed in Section 2.1; but there is no one-to-one correspondence between the two 
lists, and this justifies our separate presentation of a list of few-body problems (of course 
any many-body problem includes few-body problems as special subcases). 

Of course the number of arbitrary "coupling constants" featured by these models can 
generally be marginally increased or decreased by rescaling the dependent or independent 
variables. 

We trust the notation below to be self-explanatory; otherwise the reader should go to 
Section H, and to some relevant appendices, where these findings are more fully explained 
and demonstrated; to facilitate such a follow-up, relevant hints are provided below. 

2.1 Few-body problems 

A linearizable one-body problem (featuring an arbitrary function, see (5.22)): 

f (t) = ip{r) f(t) Af (t) . (2.1) 
Special cases 

i). The choice ip{r) = k yields 



r[t) = k r{t) A f{t) , 



(2.2) 
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which describes the motion in ordinary space of a particle acted upon by a force pro- 
portional to its angular momentum (of course the arbitrary coupling constant k could 
be rescaled away). This is a solvable equation: its explicit solution in terms of parabolic 
cylinder functions is given in Appendix A (see formulas (A. 19), (A.28) and (A. 36)). See 
(5.23). 

ii). The choice (f{r) = k/r^ yields 



fit) = k f{t)Af{t) /r^ , 



(2.3) 



which is solvable in terms of hypergeometric functions (see Appendix A, formulas (A. 19), 
(A.28) and (A.38)). See (5.25). 

iii). The choice tp{r) = k/r^ yields 



f'it) = k r{t) A f{t) jr^ , 



(2.4) 



which is the "Newton / Lorentz" equation of an electrical charge in the magnetic (radial) 
field of a magnetic monopole (or of a magnetic monopole in the electrical field of an electri- 
cal charge) . This equation is partially investigated in Appendix A, a deeper investigation 
is postponed to a subsequent paper. 

A linearizable one-body problem (featuring an arbitrary even function, and 3 arbitrary 
constants) : 



r = 2ar + br-\- c 



2f[f-f 



r \ r ■ r 



jv^ -\- r A r / (r) 



(2.5) 



See (5.17). Subcases of this model (with a = b 
subcases are solvable, including the case / (r) = 
problem (with 3 arbitrary coupling constants): 



2af+ bf+ c 



2f ■ — r (j^ ■ /r 



- c = 0) were considered above; other 
0, which yields the solvable one-body 



(2.6) 



The generic solution of this equation of motion is completely periodic, with period T 
27r/u}, if a = and 6 (c - 1) = > g. Sec (5.1), (5.2). 

Another interesting case obtains setting c = 0, / (r) = C in (2.5), getting thereby 



r = 2ar-\-br-\-CrAr 



(2.7) 



See (5.18). This one-body problem is solvable if 6 = (8/9) or 6 = 0; see Appendix B. 

An integrable and linearizable (clearly Hamiltonian) one-body problem (with 2 arbitrary 
constants): 



r = af — br'^f . 



(2.8) 



See (5.16). 

An integrable (scalar & 3-vector) one-body problem (featuring 2 arbitrary constants, 
and a constant three- vector): 



P 



2 c' 



p (p2 -Sr^) +7P- (c-f 



(2.9a) 
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r = 2c^ 



(-r^ + 3/) +7r + pC 



(2.9b) 



Of course, only for C = is rotation-invariance preserved. See (5.15). 

A solvable (scalar &; 3- vector) one-body problem (with 4 arbitrary coupling constants): 

p = a + Pp + j [p + c {p'^ -r^)] -c 3pp-3 (^f-f^ +cp {p^ -3r'^) , (2.10a) 



r = /3r + 7 



r + 2cpf 



Spf+Spf— fAf+cr (Sp^ — r^) 



(2.10b) 



See (5.12). 

A solvable (translation-invariant) two-body problem (with two arbitrary constants): 



■4±) 



r =-^Q;lr -|-r j^tG 



I / _ / 1 /\ ( _ 1 W2 . 



(2.11) 



See (5.24). 

A solvable (translation- invariant) 2-body problem (with 4 arbitrary coupling constants): 

'r^^ = [{a/2)±a\r^^ -\-[{a/2)Ta\r'~^ ±\^br-\-c 2f (j" ■ 'Fj - f ■ /r^} , 

(2.12a) 

r{t) = r<+) {t) - r<-) (t) . (2.12b) 
See (5.5). 

A linearizable two-body problem (with 6 arbitrary constants): 



r^'^ = 2 (ar''^ - c,f^'>) +Pf<'^ -^r^^^ 



+ c 



r — a r 



- I -(1) A I ^2) A 

c r- A r -I- -n ' A 



^1 \ 



-t- c r^"^ A r r^"'' A r 



r^"^ A r - ■P"^ A r 



')■ 



(2.13a) 



(2.13b) 



If a = /3 = 0, a = 3a;/2, /3 = 2u;2 orQ; = /3 = /3 = 0, 5; = w/2, with u; an arbitrary 
(real, nonvanishing) constant, this model is solvable and its generic solution is completely 
periodic with period T = 2tx j \uj\. See (5.20) and Appendix B. 

A solvable 3-body problem (with 4 arbitrary coupling constants): 



7^^'^ = (2a + cA) ^^'^ + 6r<^') 

fe=l,2,3,mod(3) 
A = f^l) . ^3) _ 

See (5.9). 



}/A , j = 1,2,3, 



(2.14a) 
(2.14b) 
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A linearizable 3-body problem (with 12 arbitrary coupling constants; note that all indices 
are defined mod(3)): 



•A 2 



{Mn+lMn+2r„+2 A Tn+i + ^ [--an,m'rm 

m=l ^ 

+ I «n+l,mMn+2fVi+2 " a„+2,mAin+l^n+l + ^ an+l,m«n+2,ferfe j Ar^n]} • (2-15) 



k=l 



See (5.27). 

A solvable 4-body problem (with 3 arbitrary coupling constants; translation-invariant if 
6 = 0): 

'f^'^=2ai^'^ + bf^^^+c Yl (-)' 

fe=l,2,3,4,mod(4) 

{r^'^ [r^'^ ■ (r-<'=+i) - r-<'=+2)) a {f<''+'^ - 7^''+'^)] } /A , (2.16a) 
A = (^2) _ ^1)^ . I^^3) _ ^i)^ ^ (^^4) _ ^1)^ _ (2.16b) 

See (5.10). 

A linearizable four-body problem (with 8 arbitrary coupling constants): 
r =ws —7s ±[ 2 \ar — ar \ + p 'n ' — p r^ ' 

+ c (r-<i) A r^'^ - r-<2) A r^'^) - c (r-<i) A r^'^ + r-<2) A r^'^) ] } /2 , (2.17a) 

•^2,±) u;2) / u;2) , ^ ^2) , 3 ^1) 

r ={75 +75 ±[2(Q;r +Q!r \ + pr^ ^ + pr^ ' 

+ c (r-<^) A r^'^ + r-<2) A r^'^) + c (r-<i) A r^'^ - r-<2) A r^'^) ] } /2 , (2.17b) 



^1) = _ ^1 ^2) ^ ^2,+) _ ^2,-) 

dX) = ^ ^;<1 -) ?(2) = ^;<2,+) , M-) 



(2.18) 



The generic solution of this model is completely periodic with period T = 2t: / \uj\\i there 
holds cither one of the two sets of restrictions on the 4 coupling constants a, a, /3, /3 
reported above (after (2.13b), and in addition there hold the two constraints 7 = 0, 7 = 
mui, with m an arbitrary (nonvanishing) integer (m 7^ 0; for m = 0, namely 7 = 7 = 0, 
(2.18) reduces to (2.13)). See (5.21). 

2.2 Many-body problems 

A linearizable N-body problem (with ZN arbitrary coupling constants): 

TV N 
ni=l ni,n2=l 

Here, and in analogous equations below, all indices are defined mod(A'^). If a„ = 6„ = 
the model is actually solvable. See (5.30). 
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A linearizable (translation-invariant) (2N)-body problem (with arbitrary coupling 
constants): 



ni=l 



.)}/2 



N 

+ E 

ni ,n2=l 

— I n ^ ' n i ' n — ' n ' n i ' n — X^n ^ ' n) I ^ ■ 



(2.20a) 
(2.20b) 



Again, if a„ = 6„ = 0, the model is actually solvable. See (5.31). 

A solvable Hamiltonian many-body problem (with 6 A'' arbitrary coupling constants): 



N 



N 



("n-ni Pni + 7n — m Qni 



■n—ni—n2—n3 



(Pni A gna) A I 



ni=l 



Tii,n2,n3=l 



AT 



(2.21a) 



Pn = - XI + 'y^-^i Pni)+ XI 



-n— m— n2— ns 



ni=l 



ni,n2,n3=l 



(2.21b) 



See (5.32). 

A scalar/vector solvable N-body problem (with 8N arbitrary coupling constants): 



N 



ni=l 
iV 

~ 3 ^ ^ -Icn— ri]^— n2 A'ni /^n2 ^''"i ' ''"s^ j ^ 

"1 ,"2 = 1 

N 

"I" ^ ^ {Cn— ni— 712— 713 Tm [/'n2 /'t13 (^"2 ' ^n3)]} 
ni,n2,n3=l 
AT 

~ ^ ^ {Cn— nj— n2— n3— n4 Cni Pn2 [Pn3 Pn4 ~ 3 (^^3 • ''714)]} j 

ni,n2,n3 ,n4=l 

AT 

''n ~ ^ ] '^/3n— m 1"ni "I" Tn— m ''ni 
ni=l 

N 

^ ^ '^Cjj_^j_jj2 3 Pjjj rjj2 + 3 Pm r^^ ^''ni A r^jj^j ^ 
ni, 712=1 

N 

"I" 2 ^ ^ {Cn— m— n2— n3 Tni /'n2 ''"3} 

»ii,n2,n3=l 
AT 

~ ^ ^ {pn—n\—n2—nz—n4 Cni ^n2 Pns A'n4 ~ (^n3 ' ^n4)]} ■ 

ni,n2,n3,n4=l 



(2.22a) 



(2.22b) 
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See (5.34). 

A solvable N-body problem (with + 3N arbitrary coupling constants): 



7n (rn ■ + CLn In A r„ 



Tn = dnrn- Oin In A f„ + {-7^ r„ A 

N 

- X] "L^"™ "L^™ + 7n r„ A [7^ {fn ■ Vra) r„ + fm + On 7n r=*n A f,„]}}} 



m=l 



/[a„(a„ + 7„ r„)] . 



(2.23) 



See (5.35). 

A scalar/vector solvable N-body problem (with AN arbitrary coupling constants and 
"nearest-neighbour" interactions) : 

Pn = (an - On+l) Cn + {ttn - ttn+l) (a^ - a^) Pn + Cn {bn Pn - K+l Pn+l) 

- {an - an) bn+1 {Pn Pn+l - Tn ' ^n+l) - 3 6n Pn - Tn ' ^n) 

+ (a„ - 2a„ + a„+i) [pn + 6„ (p^ - r^)] 

+ 6„+i {pn pn+l - Tn ■ Tn+lj + bn bn+1 [pn+1 {Pn " '^n) ~ ^ Pn fn ■ fn+l] 

-blpn{pl-3rl), (2.24a) 

fn = (an - ffln+l) (^n - ^n) r„ + C„ (6„ f„ - fn+l) 

- {an - an) bn+1 {Pn fn+l + Pn+l fn + fn/\ fn+l) 

- bn (3 Pnfn + ^ Pn fn-fn^ ^n) + (a„ - 2 a„ + 0„_i) (fn + 2bnPn fn^ 
+ bn+1 (pn+1 fn + Pnfn + fn/^ fn+lj 

+bn bn+1 [2 Pn+l Pn fn + (Pn " ^n) ^n+1 + '^Pnfn/\ fn+l] - b^ (3 p^ - r^) fn ■ 

(2.24b) 

Sec (5.36). 

A linearizable N-body problem (featuring 2A''^ + N arbitrarily assigned functions of 
time): 



2 fn. A f„ - f„ A r 



N 

Tn = ^ ^ {2 Onni (0 '"ni + ^nni {t) fni + Cni (i) 
ni=l 

-Cni (t) (f„ Afni)} 

AT 

+ X {2a„ni (*) C'na (*) (^m A fnj - C^^ (t) C„2 (t) [(f„ A f„J A r;^]} . 
ni,n2=l 

(2.25) 

See (5.38). 

An integrable N-body problem (with "nearest-neighbour" interactions): 



2 T^n ( Tr), ' r^ 



n \ ' n ' ' n ] fn \ ' fn ] l^n 
+ 7 {fn+l - [2 

fn {fn ' fn—l) fn-1 {fn-fn)]/rn-l} ■ 



(2.26) 
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See (5.39). 

A scalar/vector integrable N-body problem (translation-invariant, with "nearest-neigh- 
bour" interactions): 



Pn = C 

rn = c 



Pn (pn+i -2pn + Pn-l) - rn ' (r„+i - 2 f „ + f„_i) , (2.27a) 
{Pn+1 — '^Pn+ Pn-l) + Pn (^ n+1 2 rji 

(2.27b) 



See (5.40). 

A linearizable scalar-vector N-body problem (with 2 arbitrary coupling constants, and 
"nearest-neighbour" interactions) : 

Pn = aPn + PnPn - r„ ■ fn + C ^Pn-lPn - " - PnPn+1 + ' ^n+l} , (2.28a) 

r^i — a + pri^n ^ Pn^n ^ 

+ C jpn-l^Vi + Pnrn-1 " ^n-l A f „ - " Pn+l^n - /\ ^n+l} , (2.28b) 



where 



and 



Pn = (pnPn + ^„ • f„ j / {pi + r^) , 

Tn = {pnTn - PnTn - A f^^ / (p^ + r^) , 



Pn+1 = PnPn+l - ?'n • , 

Tn+l = Pn^n+l + Pn+l^n - A f„+i 



(2.28c) 
(2.28d) 

(2.28e) 
(2.28f) 



See (5.43). 

A linearizable N-body problem (with SAT arbitrary coupling constants, and "nearest- 
neighbour" interactions): 

Pn = (dn+l - an- K) Pn + PnPn - " + K+1 {pn+\Pn - " 

+ c„+ip„ - c„ - • Tn-lj , (2.29a) 

r„ = (a„+i - an - 6„) f„ -|- p„f„ -|- p^fn - A r„ 

+ K+l (^Pn+lTn + - Tn+l A 

+ Cn+ir„ - Cn (^Pnfn-1 + Pn-l^n + A f„_i^ , (2.29b) 

where p„,f„ are given again by (2.28c), (2.28d) and 

Pn-l = PnPn-l " ' , (2.29c) 

fn-1 = PnTn-l + Pn-l^n - Tn f\ r n-l , (2.29d) 
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Pn = {pnPn + " 4) / + ' ^n) , (2.29e) 

Tn = [pnTn - PnTn + A / (^pl + fn ' ^n) , (2.29f) 

See (5.44). 

For linearizable equations of motion involving the AT + 5 three- vectors r„ (t), f {t), g {t), 
h{t), v{t), y{t) and the 2 scalars r/(t), 6'(t), and featuring 2 AT^ + 20 AT + 42 arbitrary 
"couphng constants," we refer the reader to (5.45). 

A scalar/vector solvable N'^-body problem (with 4A^^ arbitrary couphng constants): 

N 



Pnm — 0"nm ~l~ ^ ^ [bmim Pnmi ~l~ Cfmm Pnmi] 
mi=l 

N 

^ ^ dmim2[Pnmi Pm^n {^nmi ' ^m^n) "l~ Pnmi Pm2n 2 (r^i^^j^ • rjn,2n)] 
mi,m2=l 

"I" ^ ^ dm\m2 Cmsm [Pnmi Pm2mz ('"nmi ' 'fm2mz)\ 



mi,m2,m3=l 
AT 

^ ^ d"mim2 dmsmi [Pnmi Pm2m3 PmAvn 

mi,m2,m3,m4=l 

pnmi {'^m2m[i ' frn4m) Pm2m3 {fnmi ' fm4m) Prrura {fnmi ' '^11121713) 

~l~ {ffimi A Tfn2rai) ' '^m4ra ] 1 (2.30a) 
N N 

mi=l mi,m2=l 
AT 

mi,m2,m3=l 
AT 

^ ^ dmzm4,dmim2 [ '"nmi {Pm2m3 Prmm {j'ra2m3 ' '^rmra)} 

mi,m2,m3,m4=l 

"I" fm2in3 {pnm\ Pnum (''nmi ' '"m4m)} "I" fmAfn {Pnmi Pm2m3 (''nmi ' ^m2m3)} ] • 

(2.30b) 

See (5.46). 

A solvable (Hamiltonian) N-body problem (with 2 N'^ + N arbitrary couphng constants): 

N 



#m + 4 A [qn/\ [qm A Pm]]) , 

m=l 

N 

Pn=^ {Cnm Qm - dmn + 4 A [pn A [g^ A Pm]]) , (2.31) 



m=l 



where 6nm — ^mn ? ^nm — ^mn 



For the corresponding Hamiltonian, see (F.8b). 
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3 Survey of solvable and/or integrable and/or linearizable 
matrix models 

In this section we survey solvable and/or integrable and/or linearizable matrix evolution 
equations. We focus throughout on square matrices, whose symbols are underlined. Most 
of the results reported in this section are probably well known. 

3.1 An explicitly solvable matrix evolution equation 

In this subsection we detail the explicit solution of the following matrix evolution equation: 
M = 2aM + bM + cMM-^M. (3.1) 

In writing (3.1) we implicitly assume the matrix M to be invertible; it is clear from the 
results reported below that this property is generally maintained by the time-evolution 
(3.1), unless the solution develops a singularity (this is possible in some cases, as can be 
easily seen from the solutions reported below). The 3 quantities a,b and c are arbitrary 
(scalar, possibly complex) constants (the factor of 2 in front of a is introduced for nota- 
tional convenience). In fact, cases in which a,b and c depend on the time t can also be 
treated, but for simplicity we refrain from doing so. 

The matrix evolution equation (3.1) can be solved in explicit form for arbitrary initial 
conditions. For c 7^ 1 the solution reads 



Hence M(i) is completely periodic with period T = 27r/a; , for arbitrary initial condi- 
tions, if 



with LO an arbitrary (nonvanishing) real constant, and m, n two arbitrary integers (positive 
or negative, at least one of them nonvanishing). These two conditions entail, if m 7^ 0, a 
single restriction on the 3 parameters a, h and c: 



However, in this case with m 7^ 0, the condition (3.3a) with (3.2b) entails that a and c 
cannot both be real. Perhaps more interesting is the case with m = 0, entailing (see (3.3a) 



M {t) = exp (a 7 1) { cosh (A t) 

+ A-^sinh(At) [m(0) [7M(0)]~^ - a] }^M(0) , 
7 = l/(l-c) , 



(3.2a) 
(3.2b) 

(3.2c) 




a 7 = i mu! 
A = inuj , 



(3.3a) 
(3.3b) 




(3.3c) 



a = . 



(3.4a) 



Then all solutions of (3.1) are periodic, with period 



T = 27r [6 (c-1)]-^/^ , 



(3.4b) 
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if this quantity is real. It is moreover clear from (3.2a) that in this real case, with a = 0, 
and b,c,T real, a necessary and sufficient condition to exclude that the matrix M (t) 
become singular (or, equivalently, noninvertiblc) at any time is that none of the eigenvalues 
of the (generally nonsymmetrical) matrix M(0) [M(0)]~^ be real. 

The behavior of the solution (3.2) in more general cases is sufficiently clear from the 
explicit formula (3.2) not to require any additional elaboration here. 

The solution of (3.1) in the special case 

c = 1 (3.5a) 

could be obtained by a limiting procedure from (3.2), but it deserves to be separately 
displayed. If 

a 7^ (3.5b) 
it reads 

M (t) = exp { [b/ (2 a)] [-t + exp (a t) a''^ sinh (a t)] } 

exp|exp(ai) a'^ sinh (at) (m(0) [M(0)]"^)| M(0) . (3.5c) 

This solution is periodic with (real) period T = 27r/a; , for arbitrary initial conditions, iff 

a = i muj , (3.6a) 
b = 2nmiJ , (3.6b) 

with u real (nonvanishing) and m, n two arbitrary integers, except for the restriction 
m 7^ 0, which entails that in this case a must be imaginary, see (3.6a). 

Again, we do not elaborate on the behavior of the solution (3.5c) in more general cases 
than the periodic one. 

Finally, in the case 

c = 1 , a = , (3.8a) 
the solution of (3.1) reads simply 

M{t) = exp {bt^/2) exp|iM(0) [M(0)]"^} M(0) . (3.8b) 



3.2 A class of linearizable matrix evolution equations 

In this subsection we present a technique to manufacture new (generally nontrivial) lin- 
earizable matrix evolution equations from known (possibly trivially solvable) linearizable 
matrix evolution equations. Some aspects of this technique might be new; we indicate 
below, and accordingly give due credit for, those aspects we know not to be new. 

In Appendix A we treat in some detail a specific, very simple, example, as indicated at 
the end of this subsection. To make that appendix self-contained, we explain there again 
the essence of the results presented in this subsection, since in that very simple context 
this can be done quite tersely. Hence the reader might find it profitable to scan Appendix 
A (at least its opening part), before delving in the treatment given below. 



Solvable and/or Integrable and/or Linearizable A''-Body Problems 



317 



Our treatment is based on the following Lemma. 

Assume that the N + 5 square matrices u„, /, g, h, v, y satisfy the iV + 5 matrix 
ODEs 

U u(^) m = 1 N- f /-(i) • 

5(°\ff«,5(^\...;/i(°),/i«,/i(^\...;t)=0, (3.9a) 

^{iLrn ' — m'^ ' —ni ' •••j^ 1) ^ i / ^ ' X'' ' ) • • ■ ! 

5(°),£W,£(^\...;/L(°UW,^(^),...;t)=0, (3.9b) 

G(uW,^a),^(2),...,^ = l,...,jV;/W,/W,/(^ 

5(°\£«,5(^),...;^(°),^«,^(^),...;t) = 0, (3.9c) 

l(MS!\2i(^\ni?,...,m = l,...,iV;/(0),/«,/(^),...; 

gj'\gj'\gj'\..,h('\h^'\h^'\..,t) = , (3.9d) 

^ ^ Y-j {iLrn 1 Him i i^rn i ■•■^f^ l)--)-^; ^ ^ ^ '< l} \ •••■> 
3=0 

g_^'),gj'\g_('\..,h('\hi'\h^'),..,t)vi^^=0, (3.9e) 
X;y(^)£,(u(:!),uW,^(^),...,„^ = l,...,Ar;/(0),/^^^ 

j=0 

g_('^\gj'\g_('\..,h'^'\h^'\h'^'\..,t) = . (3.9f) 

Here Un^ = dPun/dP ^ j = 0, 1, 2, with analogous formulas for f^^\ g^^\ h^^\ vj-^\ 
and the functions ^„ , F, G, F , F^-, are "scalar /matrix functions of matrices," 
namely functions built only out of their matrix arguments (whose ordering is of course 
important) and of scalars (including, possibly, the independent variable t), so that they 
satisfy the identity 

H:-^I(n(^),nW,n(^),...,m = l,...,Ar;/(0),/(i),/^^^ 
g^'\g^^\^'\..,h^'\h^^\h^^\..,t)W = 

l(Ii^"^W^^^W^,li:~^2iite,li:"^2i^^^W^, 1, AT; 

i5/-l/(0)]^,]^-l/(l)]^, ]^-l/(2)vF^ .... 

(3.10) 

where Z denotes any one of these functions. The notation employed to write these equa- 
tions is meant to suggest (in a sense which will become clear below) that the equation 
Z_ = plays the primary role in determining the matrix z , with the obvious correspon- 
dence (if Z_ = U_j^, then = u„; if = then z = f; and so on). Note moreover the 
qualitative difference among the equations satisfied by the N + 3 matrices u^, / , g , h, 
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(3.9a,b,c,d), and those satisfied by v and y , (3.9e,f) (in this latter equations the parame- 
ters Ji and J2 are of course two nonnegative integers). The different role of the N matrices 
n„, and of the 3 matrices / , g , h, will become clear below. 

Consider now square matrices C/„ , F_, G, H_, Y_, related to the N+b square 
matrices u„ , / , 5 , /i, w , y as follows: 



-1 



h = WHW 



-1 



Un = WU^W-\ l = WFW 
Un = Wr^UnK, F = W-^IW, G = W-^gW, H = W-^hW, 
V = WV, y = Y_Wr^ , 
\^ = Wr^v, Y = yW, 



with 



W = gv + lW + WyhW , 
W = WK , 



where 



K = F + GV + Y H . 



(3.11a) 
(3.11b) 
(3.12a) 
(3.12b) 

(3.13a) 
(3.13b) 

(3.14) 



Note that (3.13b) with (3.14) follows from (3.13a) via (3.11a) and (3.12a). 

To complete the definition of the (invertible!) matrix W_ = W_{t) one must supplement 
the matrix ODE (3.13) with an "initial condition," say W^(0) = W_q. It will remain our 
privilege to make an appropriate choice for the (invertible!) matrix W_q; generally it will 
be convenient to set simply W_q = 1 . 

We can now formulate the first of the two Propositions contained in our Lemma. 



Proposition 1. The N + 5 square matrices C/„ , G, F, y can be obtained from 

the N + 5 square matrices u^^ , f , g , h, v, y by algebraic operations (inversion and mul- 
tiplication of matrices ) and by solving a linear second-order matrix ODE ( which in some 
cases reduces to a first-order matrix ODE, see below). If the matrices f , g , h, v, y are 
time-dependent, this linear ODE is generally nonautonomous. 

Let us then define matrices U^^ , F^^ , G}j^ , H}j^ and Vj-^'^, ^^^Y , via the following 
recursive formulas: 



Wz = Z, ^J+^^Z = ^^Z - ^^^ZK , 
so that 

'z,K 



Z,K 



Z}^^ = Z-[Z,K], z}'^^=Z-2 
Z^^} = Z + KZ, = 1 + 2^1 KZ -\-K^Z, 

^^^Z = Z-ZK, ^'^^Z = Z-2ZK-ZK + ZK'^ 



+ [\z,k].k: 



(3.15a) 

(3.15b) 
(3.15c) 

(3.16a) 

(3.16b) 
(3.16c) 



Here and throughout \A, B\ denotes the commutator of the two matrices A and [A, B\ = 
AB-BA . 

We can then formulate the second Proposition of our Lemma. 
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Proposition 2. If the N + 5 square matrices u^, f , g, h, v, y satisfy theN + b matrix 
ODEs (3.9), then the N+h square matrices , G, I^, V_, satisfy the N+5 matrix 
ODEs 



U C/W ^7^1 m = l Ar-F[°l F^l • 

, gW , , , fl;!^] , ^[2] , = 0, 

gM, gW, GPI, ...;^[0],^[1],^[2], = 0. 

, gW , gPI , ^[01 , , ^Pl , t) = 0, 
i^m[o] ^[2] jY.^[o] ^[1] ^[2] . 

G[°1,gW,gP1,...;^[o],^W,^P]^___.^)^0^ 

, gW , gP] , ^[0] , , ^[2] , yO} = 0, 

j=0 

G[oi,GW,GP],...;^[o],fl;W,fl;P],...;t)=o, 



(3.17a) 
(3.17b) 
(3.17c) 
(3.17d) 



(3.17e) 



(3.17f) 



obtained from (3.9) by replacing, wherever they appear, the matrices Un\ f^^\ g^^\ h^-'\ 
vU)^ yU) respectively with the matrices u}i\ F^l G^l H^l V^^^, {j>r. 

Proof of the Lemma. It follows easily from the relations (3.11) and (3.12), taking into 
account (3.13) and (3.14). Indeed the matrix Riccati equation (3.13a) is linearizable via 
the position 



W = - (yh) ^ MM-^ , 
which yields for M the linear ODE 

K= (^yh + yh + yhlj {yh)'^ M-yhgvM . 



(3.18) 



(3.19) 



Hence, if the 5 matrices / , g , h, v, y are known, the matrix W_ can be evaluated, 
via (3.18), by solving the linear (generally nonautonomous) matrix ODE (3.19). Then 
the + 5 matrices L[„ , E.i G, H_, V_, Y_ can be obtained from the, assumedly known, 
A'' + 5 matrices u„ , /, g, h, v, y via (3.11b) and (3.12b). 

The first part of the Lemma is thereby proven. 

The second part of the Lemma is implied by the formulas 
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where z_ stands for n^, f,goih, and correspondingly Z_ stands for U_^, G or H_, and 

^(j) = w_ y{^> , (3.20b) 
yU) = {j}y 1^-1 ^ (3_20c) 

which, via (3.10), (3.11) and (3.12), clearly entail the equivalence of (3.9) and (3.17). 

As for the validity of these formulas, (3.20), they are an immediate consequence of 
the definitions (3.15) and of the formulas (3.11) and (3.12), as can be easily verified by 
time-differentiating them and using (3.13) and (3.14). 

The Lemma is thereby proven. ■ 

Note that we are implicitly assuming that the 2 matrices W_ and yh arc invcrtiblc, 
see (3.11), (3.12) and (3.18). A breakdown of one of these two conditions at some time 
tc will generally show up as a singularity of the corresponding solutions of (3.17). But 
note that the second condition (invertibility of yh) is invalid but irrelevant in the special 
subcase in which y h vanishes identically, see the Remark (iii) below. In this special case 
the (first-order) matrix evolution equation (3.13a) is already linear, hence there is no need 
to introduce via (3.18) the matrix M. 

Now we are ready to formulate and prove the main result of this subsection. 

Theorem. If the matrix ODEs (3.9) are linearizable, the matrix ODEs (3.17) are also 
linearizable. 

Proof. The second Proposition of the Lemma entails that the solutions of (3.17) can be 
obtained by first solving the ODEs (3.9) to determine the N + ^ matrices u„, /, h, v, y 
and by then obtaining from these, according to the first Proposition of the Lemma, the 
+ 5 matrices U^, F, G, H, V_, Y. U 



Remarks 

(i) Note that, even if the equations (3.9) satisfied by the A + 5 matrices /, g, h, v, y 
are linear (and possibly quite trivial, see examples below), the equations (3.17) satisfied 
by the A + 5 matrices F_, G, H_, V_, Y_ are generally nonlinear. Also note that, in the 
context of the initial-value problem, the convenient choice W n = i, see above (paragraph 
after (3.14)), entails that the initial conditions for the A'-|-5 matrices u^, f, g, h, v, y can 
be explicitly obtained from the initial conditions for the A"-|-5 matrices U_^, F, G, H_, V_, F 
via (3.20) with (3.15,16). 

(ii) Not all the equations (3.9a,b,c,d,e,f) need be differential, for instance (3.9b) might 
read 

L = L{^m,m = l, N;g^^^;hS^);j = 0, 1, 2, ...) , (3.21) 

with / a given "scalar /matrix function" of its matrix arguments, yielding thereby an ex- 
plicit definition of the matrix / in terms of the matrices g, h and their time-derivatives, 
and likewise an explicit definition of the matrix _F in term of the matrices L[„, G, II and 
their time-derivatives, sec (3.7b) and (3.20a); and so on (see examples below). But this 
applies only to the equations (3.9a,b,c,d), and likewise to their counterparts (3.17a,b,c,d), 
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not to (3.9e,f) and (3.17e,f) (except in the trivial cases discussed in the following Remark 



(iii) The linear homogeneous equations satisfied by v and y, see (3.9e,f), and likewise 
the (generally nonlinear) equations satisfied by y_ and Y_ , see (3.17e,f), clearly admit the 
trivial solutions v=y_ = Oory = Y_ = (the vanishing of v entails the vanishing of V_ , 
and viceversa, and likewise for y and Y_ , see (3.12)). This corresponds to reduced versions 
of the above Lemma, which are obtained by setting to zero one of the pairs of functions 
V, V_ and y, Y_ (or possibly both pairs), as indeed entailed by the special cases of (3.9e) 
and correspondingly of (3.17e), respectively of (3.9f) and correspondingly of (3.17f), with 
Y-O — 1' = for j > , respectively Y_q = 1, = for 3 > 0. In the first case, 
V = V_ = 0, the matrices g and G can also be ignored, see (3.14), and one can then 
ignore the equations (3.9e) and (3.17e) as well as (3.9c) and (3.17c). In the second case, 
y = Y_ = 0, h and H_ can likewise be ignored, see (3.14), as well as the equations (3.9f), 
t3.9b) and (3.9d), (3.17d). If both conditions hold, v = V = and y = Y = 0, then the 
only matrices that play a significant role are u„ and / , and correspondingly L[„ and F, 
and the only relevant equations are (3.9a,b) and (3.17a,b); in this special case K_ = F_, see 



We now exhibit two classes of linearizable second-order matrix ODEs satisfied by the 
+ 5 matrices C/^, F, G, H_, V_, Y_, obtained from this Theorem by making specific 
choices for the functions f/„, F, G, H_, V_, y, see (3.9) and (3.17). Once and for all, 
let us emphasize that these are representative examples, selected to display the type of 
linearizable matrix ODEs encompassed by our approach, and chosen with an eye to the 
manufacture of linearizable many-body problems, see the following Sections |^ and ^; it will 
of course be trivially easy for the diligent reader to manufacture additional examples. 

The first class of examples we consider is characterized by linearizable (in fact, linear) 
matrix ODEs of type (3.9) which read as follows: 



(iii) ). 



(3.14). 



Examples 




g + a^^^^ h + b^^^^ h 



(3.22a) 



TV 



m=l 



+ a(-^f ) g + ) g + a^-^'^) h + b'-f^^ h 



(3.22b) 



N 



m=l 



+ a(39) g + b^9g) g + a^gh) + ^{gh) ^ 



(3.22c) 



N 



m=l 



(3.22d) 
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N 



m=l 

+ ^^''^^ / ^ a^'''') ^^^'^^ /i^i^ 

TV 

+ a^''''^h^ + b^''''^h^}v, (3.22e) 

N 
m=l 

+ ^2 ^ 4% 

N 

+ y { E K'"^ + ^^'"^ + / + b^"'^ t + f 

m=l 

These equations contain N + 5 quantities of type /U, and (2A^^ + 20 A'" + 42) quantities of 
type a, b (variously decorated with lower indices and upper symbols of identification); our 
treatment would apply even if all these quantities were (arbitrarily!) given functions of 
time, but for simplicity we assume hereafter that they are (arbitrarily!) given constants. 
Then the + 3 matrices u^, f, g, h can be obtained from (3.22a,b,c,d) via purely alge- 
braic operations, and the 2 matrices v and y can subsequently be obtained by solving the 
linear nonautonomous matrix ODEs (3.22e,f). 

The corresponding (generally nonlinear!) linearizahle matrix ODEs satisfied by the 
N + ^ matrices F, G, H^, F, Y_ read as follows: 



N 



[tin,. 



K\ + 



= 2 \^P,K 

N 

+ ' 
m=l 



(3.23a) 



Un,K\-[[Un,m,m 

m=l 

+ a^r^^ {I - [F, K] } + F + ai"^) {G-[G,K]] 
h^r^ G + ai"'^) - [H, ^] } + bi"'^) H , 
^ - [[F,K],K\ 

E K^"^ [i^m - [Um.m} + Um] + ^^^^^ {Z - [Z, K]} + F 

rn=l 

+ a(-^») {g - [G, K] } + G + a^f^^ [h-\H,K]^+ b^f"^ H , (3.23b) 

/.(^)g = 2 [g,^] + [g,k]-[\g,k],k] 

N 



m=l 
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+ a(39) |g - [G, ^] } + b^39) G + a(3^) - [H, K] } + b^^^'^ H , (3.23c) 

m,K],K] 



N 

+ ' 



+ 



H,K 



+ 



m=l 

AT 



r?i=l 



f2 
-m 



+ o^''-^) {f - [F, } V 1^"^^ + 5(^5) 1^ -[G,K]y + ) 

N p 2 



m=l 



+ [£;^^]|%6K)f2^ flip's) |g- 



(3.23e) 



^iy)Y = 2YK + Yk-YK^ 

N 



m=l 

+ a^^-^) {Z -[F,K]y + h^yf^ FJ + a^ys) |g - [G, K] } V b^^^) ^2 



m=l 



+ a 



(yf) 



[t - [F, K]^ + b^y^^ F^ + a^ys) |g - [G, K] } + h^vs'^ 
+ a'^y^^ - [H,K]y + b^y^^ } . 
In these formulas K_ is of course defined by (3.14). 



(3.23f) 



The second class of examples we consider is characterized by (much simpler!) linearizable 
(or linear, see below) matrix ODEs of type (3.9) which read as follows: 



Un = u (u„, m = 1, A^; t) 
f = f{Um^iLm, m = l,...,N;t) . 



(3.24a) 
(3.24b) 



Here is a given scalar/matrix function, such that (3.24a) is a solvable or lineariz- 
able system of ODEs (see for instance the preceding Section |3.1|, or below), and / is an 



324 



M. Bruschi and F. Calogero 



(arbitrarily) given scalar /matrix function. The corresponding ODEs of type (3.17) then 
read 



+ 



[\Ur„F], F]+u 



F = f (U„„ Um - [LLm, E], m=l, N; t) . 



(3.25a) 
(3.25b) 



Note that wc arc here in the special case mentioned in the Remark (iii), see above, corre- 
sponding to the choice v = V_=0 and y = ^ = 0. 

Also note that, if / depends on the matrices (besides the matrices u^), see (3.24b), 
hence F depends on the matrices (besides the matrices C/^ ; see (3.20a) and (3.16a), 
then the ODE (3.25a) contains the second time-derivative of the matrices in the right- 
hand-side as well (from the second term). 

For the special choice 



TV 



Un=^[2 anm {t) + bnm (t) 



(3.26) 



m=l 



which entails of course that the ODEs (3.24a) are linear, the ODEs (3.25a) read 

tin = 2 [Un, f] + [Hn, l] ' [Un^ Z] , E] 

N 

+ ^ [2a„^ (t) [U^ - [[/„, F]] + bnm it) Um] ■ (3-27) 



m=l 



The factor 2 in the right-hand-sides of the last two equations is introduced for notational 
convenience, see below. In the last equation, (3.27), the matrix F_ is of course always given 
by (3.25b), with an arbitrarily chosen /. For instance, if we make for this function the 
simple choice 



N 



I {Um, ikm^ m = l, N;t)=^ [Cm {t) U^] , 



(3.28a) 



m=l 



which entails 

N 



F=^[Cm{t) Um], 



(3.28b) 



m=l 



then the linearizable matrix ODEs satisfied by the matrices C/„ read as follows: 



N 



= { 2 a„„ (t) U^ + bnm (t) [/, 



rn=l 



+ 2 Cm (t) 



LLn, Um 



+ Cm (t) 



Un, Um 



+ Cmit)[Un,Um]} 
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N 

- ^ {2anmi (*) Cm2 (*) [Hmv ^^ma] + ^mi (*) (*) [ [Hn^ U^,] , Um2] }• 

mi,m2=l 

(3.29) 

These equations deserve further elaboration, which for simplicity is hereafter restricted 
to the case with time-independent constants Gnm, bnm and Cm- Then the above ODEs 
becomes 

N 

tjin ~ ^ ^ ^ ^ 0,nm LLm ^nm Um "I" 2 C^t^ , U,^ + Cm U-n i j ^ 



N 



{2(inmiCrn2 [f/^mi 5^^7112] + Cm2 [ [^^m ^^mj i ^^^12]} ' (3.30) 

mi,m2=l 



and the corresponding ODEs satisfied by the matrices u„ read 



TV 



^ ^ [2 d^j^ itj^ + bnm ]hn\ 1 



(3.31) 



m=l 



and can therefore be solved by purely algebraic operations. Particularly simple is the 

"diagonal case" characterized by the restrictions a„m = ^nm o-n j ^nm = ^nm bn, which 
entail that the ODEs (3.31) uncouple and their solution reads 

M„ {t) = A-^exp {an t) [A„ u„ (0) cosh (A„ t) + u„ (0) sinh ( A„ t)] , (3.32a) 

An = (a^ + 6n)'^'- (3.32b) 

Even simpler is the case with a„ = a, bn = b, hence A„ = A = {a? + 6) . Then the 
evolution of the matrix W;(t), see (3.13a), (3.24b) and (3.28a), (3.32), reads simply 



W(t) = e^p{at) cosh (At) + 5 A~^ sinh (At)] W{t) , 



(3.33) 



AT AT 

with A = [cm Mm (0)]) ^ — ^ [cm Um (0)] two Constant matrices. But let us em- 

m=l m=l 

phasize that the simplicity of this case, see (3.33), has a rather trivial origin: indeed the 
nonlinearity of the matrix ODEs (3.30) in this special case a„m = ^nm d, bnm = ^nm b is 
in a way marginal, since in this case the ODEs (3.30), which can of course be rewritten as 
follows. 



Un = '2aUn+hUn + 2 



U F 



+ 



Rn, F -2aUn,E] 



[Uln, F],F], 



N 



F=^[CmUm], 



(3.34a) 
(3.34b) 



m=l 



are easily seen (by multiplying (3.34a) by c„ and summing over n using (3.34b) to imply 
a single (uncoupled) ODE for the quantity F, 



F = 2aF + bF + 



F, F 



(3.34c) 
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Note that, once this equation has been solved for F_, the ODEs (3.34a) for the matrices 
C/„ become linear (albeit nonautonomous) . 

The matrix ODE (3.34) coincides (up to a trivial rescaling) with the special case of the 
system (3.30) corresponding to = 1. Let us conclude this subsection by focusing on 
this interesting case. There is then a single matrix U_i = U_ (t) which satisfies the matrix 
evolution equation 



U = 2aU + bU + c 
or, more generally, 

U = 2 ail + bU + 



u, u 



(3.35) 



U, fim , (3.36) 



if we retain the freedom to make an arbitrary choice (rather than the special choice (3.28)) 
for the scalar/matrix function / , sec (3.24b) and (3.25b). Note that, independently of 
this choice, the corresponding evolution equation for u (t) reads simply 

u = 2au + bu, (3.37a) 

hence its solution reads (see (3.32)) 

n(t) = A"^exp(ai) [Au(0) cosh (At) + u (0) sinh(At)] , (3.37b) 
A={a^ + b)^^\ (3.37c) 

and the corresponding evolution equation for W_ coincides of course with (3.33), now with 
A = cu{0) , B_ = cu{0) (two constant matrices). 

In Appendix A we discuss in some detail the matrix evolution equation (3.36) with 
a = b = and for 2 x 2 matrices. In this special case, and for special choices of the function 
/ (including / (u) = cu which yields (3.35), the corresponding evolution equation for W_ 
can be explicitly solved in terms of known special functions (sec Appendix A). 

In Appendix B an additional analysis of the matrix evolution equation (3.35) with a, b 
and c three arbitrary (possibly complex) constants is given, and in particular it is shown 
that, if a = (3/2) iu, b = 2u)'^ , or if a = iuj/2, 6 = 0, with lo an arbitrary (nonvanishing) 
real constant, this matrix ODE is solvable and all its solutions are completely periodic, 
with period T = 27r/a;. 

The matrix ODE (3.35) is not new: its linearizability (essentially by the same approach 
as explained above) was already demonstrated by I. Z. Golubchik and V. V. Sokolov [3] 
(see also the reference in their paper to a private communication by E. V. Ferapontov 
relevant to the special case treated in Appendix A). 

Let us end this subsection by pointing out that the examples given above are focussed 
on second-order ODEs for the + 5 matrices F_, G, H_, V_, Y_ (or for a subset of 
them), obtained starting from second-or(ier gibe's for the A^+5 matrices u^^, /, g, h, v, y 
(or for a subset of them). But it is also of interest to consider first-order matrix ODEs, 
both because of their possible applicative relevance, and because such equations may also 
be connected to our main goal in this paper, namely to construct models of many-body 
problems in ordinary (three-dimensional) space. Indeed there are two ways in which such 
a goal may be realized by starting from first-order matrix ODEs of the kind yielded by the 
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technique introduced in this subsection: such equations may be appropriately interpreted 
as Hamiltonian equations of motion (rather than Newtonian equations) for a many-body 
problem, or they may be used to obtain (by time-differentiation and appropriate substi- 
tutions) new second-order matrix ODEs, which can then be appropriately interpreted as 
Newtonian equations of motion for a many-body problem. Each of these two possibilities 
is illustrated by an example, whose presentation is relegated to the two Appendices F and 
G. 

3.3 Some integrable matrix evolution equations related to the non- 
abelian Toda model 

In this subsection we consider various coupled systems of matrix evolution equations, 
which coincide, or are closely related, with the integrable evolution equations of the so- 
called nonabelian Toda lattice [4] 

Q.n = G„ [G„]"^ G^ + l [ G„+i - G„ ^ G„| . (3.38) 

Here G„ = (t) is a time-dependent square matrix, labeled by the index n, and 7is a 
("coupling") constant (possibly complex, see below), which could be eliminated via the 
scale transformation 5^ 7" or, if it is positive, via the time rcscaling t 7^1/2 t. 

In this paper we do not discuss the actual solution of this matrix evolution equa- 
tion: it suffices at this stage for us to know that this is an integrable equation [4]. 
But let us emphasize that this evolution equation should be completed by prescriptions 
"at the n-boundaries," such cts , Scty, Goit) = G^+i(t) = ("free ends") or Gq (t) = 
G^it), Gi{t) = G^+i(i) ("periodic"), if (3.38) is to hold for n = l,...,iV; of course 
these prescriptions are relevant to determine the solution. 

Another, perhaps more interesting, integrable matrix model reads 

Q =Q \Q^,-Q]-\Q -Q ,] Q ■ (3.39a) 

Here of course = (t) is again a square matrix. 

The integrability of this matrix evolution equation is entailed by its relation to the 
integrable equation (3.38). Indeed it is clear that, via the positions [4] 

^nW = [G„(t)]-'G„+i(i) , (3.40a) 
S„(t) = [G„(i)]-^G„(t) , (3.40b) 

(3.38) can be rewritten as [4] 

An =iin^n+l ~ A.n i (3.41a) 
S„ = 7 [^n-^n-l] • (3-41b) 

We now set 

n 

m=mo 
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entailing 



Q 



■n — 1 



hence, via (3.41b), 

7An = Q„ ■ 



(3.42b) 



(3.42c) 



Insertion of (3.42c,b) in (3.42c,b) yields precisely (3.39a). 

Of course also this equation must eventually be completed with "end-point" prescrip- 
tions, namely appropriate definitions for (t) and Qj^_^^ {t) (assuming again that the 
index n in (3.34a) runs from 1 to N); consistently with this ambiguity we left undefined 
the lower limit mo of the sum in the right-hand-side of (3.42a). 

Note that (3.39a) is invariant under the translation (t) — > (t) -|- Q^, = 0. 

An arbitrary ("coupling") constant c can be reintroduced in (3.39a) via the rescaling 
Q — > c Q , so that it read 



Q 



Q 



n+l 



Q 



Q -Q 



■n — 1 



Q I • 

— n J 



(3.39b) 



Note that the ansatz (t) = n U_ {t) is compatible with this evolution equation, and 
it yields for U_ (t) the simplest one of the evolution equations discussed in the preceding 
subsection, see (3.35) (with a = b = 0) and Appendix A. 

Finally, we report another integrable first-order matrix evolution equation, which is the 
simplest instance of a class of integrable nonlinear matrix evolution equations [5]. It reads 
(see (5.1) of [5]) 

iiji = C {An-iAn ~ AnAn+l} > (3.43a) 

and it yields, via the position 

An = Un [Un]'' , ILn = An LLn ^ (3.43b) 

the second-order linearizable matrix evolution equation 

Un = U,, + c {U^_, [C/„_i] ' - E„ L[„+i [/„} . 

(3.43c) 

Here of course C/„ = C/„ (t). This equation is linearizable because, to solve it, one must 
solve firstly the integrable equation ( 3.43a| ) and then a linear nonautonomous matrix ODE 
(see the second of the (|3.43bD ). 

The structure of this evolution equations entails that, if C/„ (t) satisfies it, then f/„ (t) = 
ILn (t) with r = [exp (at) — 1] /a satisfies the more general evolution equation 



(3.43d) 



This suggests that, if a = ztzcj, cj > , the generic solution of this latter equation is 
completely periodic with period T = 2-it/uj . 









1 i 










--atn+tn 






Un 


























+ C 1 tn^l 


Un^ 


^1 


tn-L 




L+1 


ILn+l 
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3.4 Some other solvable and/or integrable and/or linearizable matrix 
evolution equations 

In this subsection we present some other solvable and/or integrable and/or linearizable 
matrix evolution equations. 

A solvable equation reads as follows: 



U = al + pU + -f l^U + clTj -c [UU + 2UU + cU:'j , (3.44) 
with a, (3, 7, c arbitrary constants. Indeed by setting 

cU = V-^V, (3.45) 
one easily sees that ( p.44| ) becomes the linear, explicitly solvable, matrix ODE 

V = caV + c/3V + -fV. (3.46) 

It is easily seen that, if the 4 constants a, /5, 7, c are real, a necessary and sufficient 
condition to guarantee that all solutions of (3.46), hence as well of (3.44), be completely 
periodic, is validity of the equalities a = 7 = 0, together with the inequality c/3 < 0. 
A linearizable system of coupled nonlinear matrix ODEs reads 

tin = tin [UnV^ + E ^™ t^^™]"' Hn , (3-47) 

m 

where the time-dependent square matrices = {t) are labeled by the index n and 
the (scalar) quantities anm are arbitrary constants (there are A^^ of them, if the indices 
n, m range from 1 to A'^). Indeed by setting 

Vn = Un Unr' , (3-48a) 
U^ = VnUn, (3.48b) 

one gets for V_^ (t) the linear, explicitly solvable, matrix ODE 

Zn = ^«nmZm- (3-49) 

m 

Hence to solve (3.47) one first solves, explicitly, this linear equation with constant 
coefficients, and then the linear nonautonomous ODE (3.48b). 
The following simple matrix evolution equation is integrable 

U = c^{2U^ + CU + UC), (3.50) 

where c is an arbitrary scalar quantity (the factor 2 could of course be rescaled away), 
and the arbitrary matrix C_ is also constant (C = 0; one could of course set C = 0, or 
C_ = C 1). Indeed it can be obtained from a special case of the integrable equation (3.39b), 
as follows. Consider the following special ("periodic") solution of (3.39b): 

Q2ni^)=Mt), Q,^^,it) = B{t) , (3.51) 
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so that the 2 matrices Ait) , 5 (t) satisfy the equations 

A = c \^Ab + bA-Aa-aA^ , B_ = c \ra + ab-bb-bb^ . 

(3.52) 

Now set 

S = A + B, D = A-B, (3.53a) 
A = {S + D)/2, B = {S-D)/2, (3.53b) 

so that 

S = ci^^AB + AB + BA + BA-AA-AA-BB-BB^ , (3.54a) 
D = c^^AB-AB-BA + BA-AA-AA + BB + BBj . (3.54b) 

It is now clear that (3.54a) can be integrated once to yield 

S = -c {D^ + 0} , (3.55a) 
while (3.54b) can be rewritten as follows: 

D = -c\^SD + DSj . (3.55b) 

But, via (3.55a) and the identification 

D{t) = U{t) , (3.55c) 

this last equation yields precisely (3.50). Note that the initial data for (3.50), U_q and 
C/g, determine, via (3.55c), i2o ^iid D.q', moreover (3.55a), with C_ assigned (arbitrarily!) 
and f/g (hence ^q, see (3.55c) given, determines S_q, while S_q can be assigned arbitrarily, 
consistently with the translation- invariance of (3.55). Prom the initial data for D(t)and 
S (t) one obtains, via (3.55b), the initial data for A (t) and B (t), namely A{0), A (0), B (0), 
B_(0); one then solves (3.52), or equivalently, via (3.51), (3.39b), and in this manner one 
finally gets A(t), B (t), hence, via (3.53a) and (3.55c), the solution U_{t) of (3.50) ( with 
c an arbitrarily assigned scalar constant and Q an arbitrarily assigned constant matrix). 
A solvable matrix evolution equation reads as follows: 

U = 2a (UU + UU) -2a^U^ -4bU^ + 3bU-2b'^U 

-2a (aU^ + bU-Uy^^ u(aU^ + bU-Uy^^ , (3.56) 

where a, b are 2 arbitrary scalar constants. 

To demonstrate the solvability of this equation, we start from the matrix evolution 
equation 

M = aM^ + bM, (3.57) 
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whose solution reads 



M(t) = |[M(0)]~Vxp(-6t) + (6/a) [exp - 1]} ^ 
Now set 



(3.58) 



(3.59) 



a position which is clearly compatible with (3.57) and (3.58), and which yields for the 2 
matrices A (t) and S (t) the equations 



A = a {A^ + B^) +bA, B = a{AB + BA) + bB, 
the first of which can be solved for S , yielding 

1/2 



B 



A/ a - A^ - bA/a 



(3.60) 



(3.61) 



The matrix ODE (3.56) is then easily obtained, via the identification ^(t) = A{t)^ by 
time-differentiating the first, and using the second, of the (3.60), as well as (3.61). 
A somewhat analogous integrable matrix ODE reads as follows: 



U = 2c icU^ + 



cU + U 



1/2 



u 



cU + U 



1/2 



(3.62) 



The starting point to get it is the integrable equation (3.43a), in the special ("periodic") 
case with A^±3 = A^. Indeed it obtains via the following steps: (i) consider the special 
case of this integrable equation characterized by the additional restriction AQ+A1+A2 = 0, 
whose compatibility with the equations of motion is easily verified; (ii) use this restriction 
to eliminate, say, A2; (iii) set Aq + A^ = S_, Aq — A^ = D; (iv) express S_ via D and D, 
thereby obtaining a second-order ODE for D that does not contain S_; (v) finally make 
the identification = —2U_. 

Next, we report the solvable matrix evolution equations 



N 

m=l 
"I" dnm Cm 



U V 



U-m' —m 



} 



(3.63a) 



where the matrices are obtained, in terms of the matrices and their time-derivatives 
U_^, by solving [6] the matrix equations 



N 



m=l 



N 



m'=l 



(3.63b) 



Here the N [3 + 2 A'") constants a„, Cn, dnm, bnm are essentially arbitrary, while the 
A^^constants dnm are the matrix elements of the matrix D, which is the inverse of the 
matrix D having matrix elements dnm- The solvability of these equations of motion is 
demonstrated in Appendix H. 
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Next, we report here two solvable and one linearizable matrix evolution equations, 
referring for their derivation (and the display of more general cases) to Appendix E. 
Solvable equations: 



U = c + aaU+{a-a) U + baU^-2bUU, 

U-n = V^n - CLn+l] C + [a„+i - a„] [Un - CLn] Hn + Cn [bnUn ~ ^n+l Hn+l] 
+ [cin — On] bn+1 lLn+1 ILn " 2 6„ [/„ L[„ 

+ [a„ + o„+i — 2a„ + bn+1 
Linearizable equation: 

Kn = K+1 -an- bn] M„ + {m„ [MJ"^ + bn+1 Kn+1 [Kn+l] Kn 



(3.64) 



(3.65) 



M 



n-l 



+ Cn+lM_^ - Cn Kn UnV Kn-1 

Finally we report the first-order integrable "Nahm equations" [7] 
Mn = c [M„+i, M„+2] , n = 1, 2, 3 mod (3) , 



(3.66) 



(3.67) 



where the constant c could of course be rescaled away (see below). A simple way to obtain 
from these equations a set of 3 coupled linearizable second-order matrix ODEs is via the 
positions 



cKn = Mn + XI '^'^."i^^m' n = 1, 2, 3 mod (3) , 



(3.68) 



m=l 

which transform (3.67) into 



m=l 



+ 



i+l^n+1 + X '^n+l,mUm^ l^n+2lLn+2 + X ^ri+2,mlLr. 



m=l 



m=l 



n = 1, 2, 3 mod (3) 



(3.69) 



Clearly we are assuming here, for simplicity, that the 3 quantities as well as the 9 
quantities a„^m, are (arbitrary) constants. These equations, (3.69), are categorized as 
linearizable, since to solve them one must first solve the integrable ODEs (3.67) for M „ (t) 
and then the linear (generally nonautonomous) ODEs (3.68) for (t). 



3.5 Duplications, quadruplications, multiplications 

In this subsection we describe some tricks whereby, from an evolution equation for one 
or more matrices, one can formally get (coupled) evolution equations involving a larger 
number of matrices. Most of these tricks can be iterated, or used sequentially, to get 
equations involving more and more matrices. These tricks are introduced and explained 
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below by showing how they work on some specific examples; they are of course applicable 
much more generally. 

Complexification. A standard method to double the number of (dependent) variables 
is to complexify. For instance setting 

U=2L + iY,a = a + ia,b = P + iP, c = j + i^, (3.70) 

one gets from (3.35) the linearizable system of 2 real coupled matrix evolution equations 



i: = 2 (aX-aYj + (^pX-pYj + 7 { 



} 



{ 



+ 



X 



Y = 2 (aX + aYj + (/3X+/3y) +7 { [k, X] + [t, 2^ } 

+7 • 



(3.71a) 
(3.71b) 



Note that there may well be other motivations to complexify than merely the purpose to 
double the number of dependent (matrix) variables under consideration; for instance we 
have seen above that there are cases with complex coupling constants when all the solutions 
of a matrix evolution equation are completely periodic (see above, second paragraph after 
(3.37c), and Appendix B). 

Association with other (solvable and/or integrable and/or linearizable) equations. For 
instance consider the solvable equation (3.1), which we rewrite here with a trivial nota- 
tional change: 



D = 2aD + bD + cDD~^D . 
Now associate to this the (trivially solvable) matrix evolution equation 

S = aS, 
whose general solution reads 

S{t) = S{0) + S{0) [exp(ai) - 1] /a . 
Then set 

Ui+) (t) = \S{t) + D (t)] /2, (t) = \S{t)-D (t)] 12 , 

S (t) = (t) + (t) , D{t) = it) - [/(-) (t) , 



(3.72) 



(3.73a) 



(3.73b) 



(3.74a) 
(3.74b) 



and thereby obtain for (t) and lA ^ (t) the system of (clearly solvable) coupled matrix 
equations 

U}^^ = [{a/2) ± a] U^^^ + [{a/2) T a] U^~^ ± l^bD + cDD'^ /2 , (3.75) 

of course with D (t) defined in terms of f/*-^-' {t) and lA^^ {t) by the second of the (3.74b). 
An important advantage of this kind of trick is to yield translation-invariant models; 
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indeed (3.75) is invariant under the translation lA"^^ (t) — >■ lA"^^ (t) + 1£q , with U_q an 
arbitrary constant matrix. 

Quadruplication. Given an evolution equation satisfied by the matrix M (i) , it is easy 
to get equations for 4 matrices, say A{t) , 5 (t) , C_{t) , D (t), via the position 

The following related formulas are worth recording, for future use: 

1 _ / {A-BD-'cy' -A-'B [D-CA-'B)-' \ 

~[-D~'C{A-BD-^C)-' {D-CA-'B)-' J' ^^'^^^ 

j^^_(AA + BC AB + BD\ 

--~\CA + DC CB + Dd) ' ^^-^^^ 
MM-'M= (If), (3.79a) 

A = (^A - BD-^Cj [A- BD-^C)~^ A 

+ [b- AA-^ Bj {D - CA-^ B)~^ C, (3.79b) 
B=(^A- BP-^ c) {A-BD-^C)'^ B 

+ (b- AA-^ B}j {D- CA-^ B)~^ D , (3.79c) 

c=(q-dd-^c^ {a-bd-^0~^ a 

+ [D-CA-^Bj{D-CA-^By^C, (3.79d) 

D=(c- DP-^ {A-BP-^C)~^ B 

+ [d- CA-^ Bj {P - CA-^ B)~^ P. (3.79e) 

The special case of these formulas with jB = Q and ^4 = ^ is also worth recording 
{^^duplicaUon"): 

M=(| j) , (3.80) 



^-,_( {A-BA-'By' -A-'B{A-BA-^By'\ 

~[-A-'B{A-BA-^By' {A-BA-^B)-' J ^^'^^ 



1 -A-^B\ ( {A-BA-^B)~'^ \ 

-A-^B 1 J 1^0 [A-BA-^Sy^ )' 



(3.82) 



MM-^M=(^^ I) , (3.83a) 

A=(A-BA-^B}j {A-BA-^B)~'^A 

+ (p_-AA-^B^{A-BA~^B)'^B (3.83b) 
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R= [a- BA-^ bJ {A-ba~^b) ^ b 

+ {b- AA~^ B}j {A - BA~^ B)~^ A- 
For instance insertion of the duplication formula 



U{t) 



into (3.1), yields the following system of 2 coupled matrix equations: 



(3.83c) 



(3.84) 



G = 



(3.85a) 

(3.85b) 
(3.85c) 
(3.85d) 



Multiplication. Finally, let us indicate two tricks whereby, from an equation involving 
one matrix, one can obtain equations involving many matrices. Again, wc illustrate each 
of these techniques via an example, but it is clear that they can be used quite generally 
(also applied taking as starting point equations involving more than a single matrix). 

First trick. Let us again consider the linearizable matrix equation (3.35) satisfied by 
the matrix [/ (t), and let us set 



(3.86a) 



as well as 



j=i j=i j=i 



(3.86b) 



where the J quantities rjj are the elements on an Abelian algebra satisfying the multipli- 
cation law 

Vj Vk = Vk Vj = r]j+k; j, k = l, J, mod (J) . (3.87) 
Here J is an arbitrary positive integer. A standard representation of this algebra is 

r]j = exp(2 7rij/J) . (3.88) 
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It is then clear that insertion of the ansatz (3.86) into (3.35) yields the following system 
of Jcoupled equations for the J matrices U_j (t) : 



■J J 



Hkl ) ll.k2 



k=l 



ki,k2=l 



(3.89) 



of course with all indices defined mod(J). 

Now assume that U (a, b, c; U{0) , U{0) ■ is the solution of (3.35) corresponding to 

the given constants a,b,c and to the initial conditions U_{0) , fZ(0); then (see Appendix 
C for a proof) the solution of (3.89) with initial condition U_j (0) , U_j (0) , j = 1, J is 
given by the (rather explicit!) formula 

J J 
If. (t) = ^{ exp i-2TTij k/J)U [ a = ^ afc^ exp(2 7riA;i k/J) , 



k=l 



fei=i 



b = bk2 exp (2 TT i A:2 /c/ J) , c = c^^ exp {2i: ik^k/ J) ; 



k2=l 



k3=l 



U{0)=Y. (0) exp (2 TT i A;4 k/J) , 

k4=l 
J 

L[(0)= 5]C/fe,(0)exp(27rzA;5A;/J)];t} . 

feB=l 



(3.90) 



Hence (3.89) is a linearizable equation, no less than (3.35) (indeed solvable for 2x2 
matrices if all the constants aj and bj vanish, aj = bj = 0; see Appendix A). 

Let us however emphasize that this trick, by its very nature, generally yields coupled 
equations that can be uncoupled by a linear transformation. Indeed it is easily seen (using 
(C.8)) that the coupled ODEs (3.89) are transformed into the uncoupled ODEs 



U„ = 2 a„ £„ + bn Hn + c, 



U-n^ —n 



via the following transformation among tilded and untilded quantities: 

J J 
Zn = J~^ Zj exp (2 TT i n j/ J) , Zn = Zj exp (— 2 tt i n j/ J) . 



(3.91) 



(3.92) 



Second trick. This is a natural extension of the quadruplication technique mentioned 
above. To illustrate it, we start from the following solvable matrix ODE: 

M = A + MB + MC + MDMC - MDM - 2MDM - MDMDM . (3.93) 

Here A, _B, C_, D are 4 constant matrices (which can be chosen arbitrarily, except for 
the condition that D be invertible, see below), and the solvability of (3.93) is demonstrated 
by introducing the matrix V_ (t) via the position 



M = D-^V-^V 



(3.94) 
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and then noticing that V_{t) satisfies the following third-order linear matrix ODE (with 
constant coefficients, hence solvable): 

V = VDA + VB + VC. (3.95) 

Note that, up to trivial notational changes, (3.93) reduces to (3.44) in the special case 
A = al, B = (31, C = 7 1, D = cl. 

We now assume the matrix M to be itself made of matrices, and we denote, with 
self-explanatory notation, its matrix elements as U_nmj likewise we denote as anmh bnmX, 
Cnm i) dnm i the matrix elements of the matrices A, B^, C_, D (the choice of multiples of the 
unit matrix is merely for simplicity). In this manner we obtain from (393) the following 
solvable system of N'^ coupled matrix ODEs for the N'^ matrices U_nm assume here 
the indices n,m to range from 1 to A^): 

N 

U-nm ~ '^nm.i+ ^ ^ ^niimLLnmi ~^ '^mimtLnmi 
mi=l 

TV 

mi,m2,m3=l 
N 

~ ^ ^ dmim2 
mi,JTi2=l 
N 

+ ^ ^ dmim2 dm^m^ ILnmi lLm2m3 ILrrum • (3.96) 

nil ,m,2 ,m3 ,m4=l 

4 Reduction of matrices to 3-vectors 

In this section we review various convenient parameterizations of matrices in terms of 
3-vectors (and, if need be, scalars). Obviously some of these results could be trivially 
extended to vectors of higher, or lower, dimensionality than 3; but, as already mentioned, 
in this paper we prefer to focus exclusively on 3-dimensional (ordinary!) space. We always 
denote (square) matrices by underlining their symbols, and 3-vectors by superimposed 
arrows. Because of the structure of the matrix evolution equations of Section ^, we are 
particularly interested in parameterizations of matrices in terms of one or more 3-vectors 
(and possibly some scalars), which belong to one (or more) of the following three categories. 

We term parameterizations of type (i), for invertihle matrices, those which are pre- 
served under the operation M M, namely are such that, if both M and M are so 
parametrized in terms of one or more 3-vectors, the combination M M admits the 
same parameterization in terms of (appropriately defined) 3-vectors. 

We term parameterizations of type (ii) those which are preserved for commutators, 
namely are such that, if both M and M_ are so parametrized in terms of one or more 



3-vectors, the commutator 



M, M 



M_M_ — M_M_ admits the same parameterization in 

terms of (appropriately defined) 3-vectors. 

Finally we term parameterizations of type (in) those which are preserved under the 
product operation, namely are such that, if both M and M_ are so parametrized in terms 
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of one or more 3-vectors, the product MM admits the same parameterization in terms 
of (appropriately defined) 3-vectors. Obviously parameterizations of type (in) are also of 
type (a), and, for invertible matrices, of type (i) as well. 

We indicate with the symbol = the one-to-one correspondence that the parameteriza- 
tion under consideration institutes among matrices and 3-vectors. For instance the more 
common parameterization we use, for (2 x 2)-matrices, reads 

M = pi -l-if- a , (4.1a) 

where p is a scalar and the 3 matrices ct^,, g_y, are the standard Pauli matrices, 

(4.2) 



1 

1 



i 
-i 



1 
-1 



(Note that in the following the unit matrix 1 = 
this case, in correspondence to (4.1a), we write 

M= {p, r) , 



Q ^ sometimes omitted). So in 



(4.1b) 



and, via standard calculations (see Appendix D, where we report for convenience a number 
of standard formulas involving a-matrices) , we also have 



M-i=(p,-r)/(p2 + ,.2) , 

M M = {^pp — r-f, pf + pr — fAf^ , 

M M={p p~p + p {r -Pj - p {r -Pj + ~p 



(4.1c) 
(4.1d) 



+ 



A r , r 



p p + [r ■ r 



Pj 



PP + 



• Pj 



+ r 



PP + [f ■ r 



+ 



pr A f — pf Af + pf Ar 



(4.1e) 



The formula (4. Id) shows that this parameterization belongs to type (in). 

Let us then restrict consideration to the class of traceless (2 x 2)-matrices, that admit 
the parameterization 



M 



ir ■ a 



(4.3a) 



(namely the special case of (4.1) with p = 0). The formulas written above entail that 
this parameterization belongs both to type (i) and to type (ii), but not to type (Hi). The 
relevant formulas read: 



M=f 



M' 



-1. 



MM-^M= 
M, m1 = - 



2 f ■ rjj — r • f j /r^ , 
2r A r . 



(4.3b) 
(4.3c) 

(4.3d) 
(4.3c) 
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The next parameterization we consider is, in terms of 2 three-vectors, for (invertible) 
antisymmetrical (4 x 4)-matrices. It reads 



M 








y{l) 




\ 


























yi2) 


_^(2) 





/ 



(4.4a) 



which we write 
It is then easy to verify that 



(4.4b) 
(4.4c) 



42) 
r 



^(2).^-<l) 



The last formula shows that this is a parameterization of type (i) (provided the 2 three- 
vectors f(^) , f(^) are not orthogonal) . 

The special case of this parameterization with r^^^ = r, r^^) = A r , 



M 



/Ox y ^ \ 

—X Az —Ay 

— y —Xz Ax 

\ — z Ay —Ax / 



(4.5a) 



with A an arbitrary (nonvanishing) constant, is also of type (i), yielding 



MM''^M= 2r (r ■ -Fj - f (r -fj jr^ 



(4.5b) 
(4.5c) 



but these two formulas merely reproduce (4.3b) and (4.3d). 

A parameterization of (3 x 3)-matrices in terms of 3-vectors takes the natural form 



M 



X(2) y(2) 

x(3) y(3) 



.(1) 

.(2) 
,(3) 



(4.6a) 



which we write as 

M = 
Then clearly 




(4.6b) 



(2) 



(3) 



,(2) „(3) 



.(2) 



u 



y 

(3) 



(4.6c) 
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which we can also write 

= u(^\ , (4.6d) 

where the 3 three-vectors u^^^ are defined, in terms of the 3 three-vectors r^^'^ , so that 

U^i) ■ji'^) =5jk\ j,k = 1,2,3 , (4.6e) 

which also entail 
3 

^^0)^0) = 1^ (4.6f) 

3 3 

^ ^0) x(^) = x^^^ = , (4.6g) 

as well as the analogous 6 equalities obtained by cyclic permutations of the x, y, z com- 
ponents of the 3- vectors n^-') = {u^^ ,'u>f'^ and f^^^ = (^x^^\y^^\ z^^^) . 

An explicit definition of the 3-vectors u^^^ reads 

^(i) ^ ^i+i) ^ j = 1^ 2, 3, mod (3) , (4.6h) 

A = r<^) • r<^) A r<^) . (4.6i) 

Note that A coincides (up to a factor 1/6 , and possibly a sign) with the volume of the 
tetrahedron of vertices 0, r^^^ , r*-^-* , r^^^ . 

Hence for this parameterization we can write 

M=(r-<^), i = l, 2, 3) , (4.61) 

and 

MM-^M= (v^^\ j = 1, 2, 3) , (4.6m) 



with 



k=l,2,3,mod{3) 



r /A , (4.6n) 



with A defined by (4.6i). 

The formula (4.6m) with (4.6n) holds a fortiori if M = M, hence it is clear that this 
parameterization is of type (i). 

A special case of this parameterization is obtained by replacing the 3 three- vectors f^^^ 
as follows: 

3 
fc=i 
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Then the treatment given above remains applicable, with the constraint on the (new) 
input vectors r to have zero sum, 

3 



(4.7b) 



It is then clear that the (new) vectors v^''\ see (4.6m,n), 



E 

fc=l,2,3,mod{3) 



r • r A r 



r /A , 



(4.7c) 



also satisfy the condition to have zero sum, 

3 



. 



(4.7d) 



Hence this parameterization is also of type (i). It has the advantage to yield, in terms of 
the original 3- vectors r^^\ translation-invariant equations. 

The next parameterization we consider is, for (4 x 4)-matrices, in terms of 4 three- 
vectors. It reads 



M 



( p 




yd) 


^(1) 


\ 


p 






^(2) 




p 


^(3) 




^(3) 




\ p 


^(4) 


y(4) 


^(4) 


/ 



(4.8a) 



which we denote as follows: 

M= (p, r<^\ J = 1,...,4 



(4.8b) 



It is a matter of standard vector and matrix algebra to obtain the corresponding formula 
for the (4 x 4)-matrix M M~^M: 



MM- 



-^M={fp-\ v(^\ j = l,...,4) , 



fe=l,2,3,4,mod(4) 



/A 



(4.8c) 

(4.8d) 
(4.8e) 



Note that the quantity A defined by this formula is translation-invariant and coincides, 
up to a factor 1/6 and possibly a sign, with the volume of the tetrahedron of vertices 

r<^), J = 1,2,3,4. 

These formulas entail that this parameterization is of type (i). 
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The last parameterization we consider is applicable to antisymmetric (3x 3)-matrices. 
It reads 




M= I -X z \ , (4.9a) 

\-y -z / 

which we denote as follows: 

M=f . (4.9b) 
The following formulas are then easy to verify: 

—y z X z 

M_M_= I —zy —xx — zz —xy | , (4.9c) 

—y X —y y — z z 




M, M = - fAf , (4.9d) 
M2'"+i=(-)"^ r2"^f, m = 0, 1, 2,... . (4.9e) 

These formulas entail that this parameterization is of type (ii) (note, however, that 
(4.9e) entails that this parameterization is also preserved for any odd power of the matrix 
M ). 



5 Results 



In this section we report, and in some cases briefly analyze, several solvable and/or inte- 
grable and/or linearizable few- and many-body problems in ordinary (3- dimensional) space, 
which are easily obtained by combining the results described in the preceding two sections. 
We do not display the calculations that underlie these findings; they are quite straightfor- 
ward, although sometimes tedious. Some of these results have been already synthetically 
highlighted in Section ||. 

The presentation is conveniently split in two parts, dealing respectively with few-body 
problems (Section |5.1|) and many-body problems (Section |5.2[) . The clarity thereby gained 
overcompensates for the minor repetitions entailed by this separation. 

Clearly the results displayed herein are far from exhausting, or systematically present- 
ing, all the solvable and/or integrable and/or linearizable N-body problems that can be 
obtained by straightforward applications of the results and methods described in the pre- 
ceding two sections. The diligent reader will have no difficulty in manufacturing many 
other solvable and/or integrable and/or linearizable A'^-body problems, after having be- 
come familiar with the methodology to obtain and investigate such problems, as described 
below (indeed, several hints in this direction are provided). 



5.1 Few-body problems 

In this subsection we consider several solvable and/ or integrable and/or linearizable few- 
body problems, obtained by applying some of the parameterizations of Section |^ to some 
of the solvable and/ or integrable and/ or linearizable matrix ODEs of Section ^. 
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A simple solvable one-body problem is obtained by applying the parameterization (4.3) 
to the solvable matrix evolution equation (3.1). Its equation of motion reads 



2 a r + bf+ c 



2 r ( r • f 



r \ r ■ r 



(5.1) 



Its general solution reads (see (3.2)) 



f{t) = f (0) L(+) {t) + iB (t)] - i f (0) ip'--^ {t) /C 



(±) 



(t) = exp (a 7 1) { [cosh {At) + {A + i C) A"^ sinh (A t) 



± [cosh {At) + {A-i C) A-^ sinh (A t)] ^ } /2, 
7 = l/(l-c) , 

A=[a' + h{l-c)]^'^ , 



A 
B 
C 



a+(l-c) f(0)-f(0) /[r(0)] 



f (0) • f(0) / f (0) Af (0) 



f(0)Af(0) / r(0) 



(5.2a) 

(5.2b) 
(5.2c) 

(5.2d) 

(5.2e) 

(5.2f) 

(5.2g) 



Here the symbol \v\ indicates the modulus of the 3- vector v, so that \v 
(irrespective of whether the 3- vector v is real or complex). 

The behavior of this system can be read from this explicit formula (see also Section 3.1). 
Of course for the model (5.1) to admit of a "physical" interpretation the constants a, 6, c 
must be rea/, as well as the initial conditions r (0) , r (0); then 93'^"'") {t) is real, (f^~^ (t) is 
imaginary (see (5.2b)), and of course f(t) is real. Note that all solutions of this problem 
are periodic, with the period given by (3.4b), if a = and 6 (c — 1) > 0. 

A solvable 2-body problem obtains from (5.1) by complexification: 



f (i) = f(^) (t) + if(^) (t) , a = a + ia, b = f3 + ij3, c = -i + i^. 
Its equations of motion read as follows: 



(5.3) 



■41) 
r 



+ 



jRW-jR(^) 

(f?(l) . f(l) _ f{2) . ^2)^2 ^ ^ ^ _^2))2 ' 



+ 



(f(l) . f(l) _ f(2) . ^2) j2 _^ 4 ^^1) . ^2)^2 ' 



where 



^(1) ^ (^l) . ^1) _ ^2) . ^2)^ ^1) ^ 2 (^f=<l) . f=<2)^ p<2) 



(5.4a) 



(5.4b) 



(5.4c) 
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^(2) ^ ^^l) . ^1) _ ^2) . ^2)^ ^2) _ 2 (^^1) . ^-<2)^ ^1) ^ (5 4^) 



f . f —J- .J- \ > + Z \r ■ r \ -n ' , (5.4ej 

1,(1) 1^2) :j;2)\ /2) o /"-(I) /I) 



j r^^^ - 2 (^r • r ' j f^^^ . (5.4f) 

An "unphysical" aspect of these equations of motion is the appearance of certain com- 
ponents of the force which are independent of the coordinate and velocity of the particle on 

^2) ~ 

which the force acts (we refer for instance to the terms —2ar and —f3r^'^^ in the right- 
hand-side of (5.4a)). This phenomenon is characteristic of several equations considered in 
this paper, and will not be highlighted again in the following. 

It is obviously easy to get from (5.2) and (5.3) the explicit solution of these equations 
of motion; again we forsake any discussion of their behavior (see Section |3.1| , and note 
that conditions on the "coupling constants" are identified there which are sufficient to 
guarantee that all solutions of this 2-body problem be completely periodic). 

A solvable translation-invariant 2-body problem is obtained by applying the reduction 
(4.3) to the coupled matrix evolution equations (3.75) (instead of (3.1)). It reads 

= [(a/2) ± a] r-^^ + [(a/2) T a] r '^ 

ib|6r + c 2f(^f-'F^—f(^f-f^ /r^l /2 , (5.5a) 

where 

f (i) = r<+) (t) - r<-) (t) . (5.5b) 

These equations of motion are translation-invariant. Their solution is given, via (5.5b) 
and 

s{t) = r<+) (t) + r<-) (t) , (5.5c) 
by (5.2) and (see (3.73b)) 

s (t) = s (0) + k{0) [exp (a t) - 1] /a , (5.6a) 

which correspond to the trivially solvable equation of motion satisfied by s{t) , 

s = as. (5.6b) 

A solvable translation-invariant 4-body problem obtains from the previous one by com- 
plexification (see (5.3)). It is left for the diligent reader to write out the corresponding 
equations of motion (a tedious task) and to analyze the behavior of their solutions (it is 
easy to identify restrictions on the "coupling constants" which guarantee that all solutions 
remain confined or are completely periodic). 
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All the models presented above can be generalized by using the parameterization (4.1) 
rather than (4.3). Every 3- vector r{t) gets then associated with a scalar p (t). The diligent 
reader will have no difficulty to obtain these results. 

Another solvable 2-hody problem is obtained by applying the parameterization (4.4) to 
the matrix evolution equation (3.1). Its equations of motion read 



r 



2ar 



_L -(1) _L /^^) 
+ b-n ' + c -i r 



r^i) \(-- 



r 



^2) 



/ |^f(l) . f(2) 



(5.7a) 



■42) 

r 



^2) 



/(r-<2).r-<i)) . 



ii2) 

r 



{f^'' . .-<1)) + {f'^ . .-<2))] _ ,-<2) p . ,^1)) } 



(5.7b) 



The task of providing an explicit solution of these equations of motion is purely alge- 
braic, but rather tedious. We leave it as an exercise for the diligent reader. 

The reduction f^^^ = r, f^^^ = A r is clearly compatible with these equations of motion, 
and it yields back (5.1). 

A solvable 4-body problem can be obtained from (5.7) via the complexification (5.3); 
the task of writing the corresponding equations of motion, and of identifying the cases 
in which all solutions are confined, multiply periodic or completely periodic, is left as an 
easy exercise for the diligent reader. 

A solvable translation-invariant 4-body problem can be obtained, as above, by applying 
the parameterization (4.4) to the coupled matrix evolution equations (3.75) (rather than to 
(3.1)), and a solvable translation-invariant 8-body problem can subsequently be obtained by 
complexification. It is an easy task (left for the diligent reader) to write the corresponding 
equations of motion, and to analyze them at least to the extent of identifying the cases in 
which all solutions are confined, multiply periodic or completely periodic. 

Next, let us consider the solvable 3-body problem that is obtained by applying the pa- 
rameterization (4.6) to the matrix evolution equation (3.1). The corresponding equations 
of motion read 



r = 2ar + b'r-'^ + c 



E 



{r [r 



fe=l,2,3,mod(3) 



/A, j = 1,2,3, 



with 

A = . ^2) ^3) _ 

Via the simple transformation 
r<^') (t) -^exp(At) r<^) (t) 
this model takes the more general form 



(5.8a) 



(5.8b) 



(5.9a) 



(2a + A (c-2)) + {b + 2Xa + X^ (c - 1)) 
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r 



k=l,2,3,mod{3) 



1,2,3, 



(5.9b) 



with A always defined by (5.8b), and it turns out, remarkably, to be closely related 
to (and to some extent more general than) an already known solvable 3-body problem 
[2]. It is an amusing task (left for the diligent reader) to compare the explicit solutions 
given in [2b] with those entailed by the treatment given herein. Also interesting are 
the various generalizations of this model that can be introduced in close analogy with 
the treatment given above (by complexification and by association with another solvable 
model so as to get translation-invariant equations); in this manner one gets solvable 6- 
and 12-hody problems, and easily identifies the cases in which their solutions are confined, 
multiply periodic or completely periodic (for instance, all solutions of the 6-body problem 
obtained from (5.9b) by complexification are completely periodic, with period T = tt/u), 
if A = 6 = 0, a = ±iu,(jj > 0). 

The equations of motion (5.9b) (as well as, a fortiori, (5.9a)) are clearly consistent with 

^ 3 . 

the restriction that the "center of mass" ^ = (1/3) ^"'^ be at rest in the origin, ^ = 

.The special solutions of (5.9b), or (5.8), that fulfill this constraint, correspond to the 
solution of the solvable translation-invariant 3-body problem that is obtained by formally 
replacing in (5.9b) or (5.8) every 3-vector f^^^ with r"''^ = r^^^ — R . 

Next we consider the solvable 4-body problem that is obtained by applying the pa- 
rameterization (4.8) with p = 1 to the matrix evolution equation (3.1) with 6 = 0. Its 
equations of motion read 

■ilj) „ liJ) sr-^ , .k 

r = 2ar + c \~) 

fe=l,2,3,4,mod(4) 

A = (f(2) - r-<^)) • (r-<3) - r^i)) A (r^^^ - r<i)) . (5.10b) 



Again, this model (which, for 6 = 0, is translation-invariant) turns out, remarkably, to be 
closely related to (and to some extent more general than) an already known solvable 4- 

body problem [2]. It is an amusing task (left for the diligent reader) to compare the explicit 
solutions given in [2b] with those entailed by the treatment given herein. Also interesting 
are the various generalizations of this model that can be introduced in close analogy with 
the treatment given above (by complexification and by association with another solvable 
model so as to get translation-invariant equations); in this manner one gets solvable 8- 
and 1 6-body problems, and easily identifies the cases in which their solutions are confined, 
multiply periodic or completely periodic (for instance, all solutions of the solvable 8-body 
problem obtained by complexification from (5.10) are completely periodic with period 
T = n/uj, if 6 = 0, a = ±iuj, lu > 0). 

Next, let us exhibit the solvable 2-body problem that is obtained by applying the sim- 
ple parameterization (4.3) to the solvable system of 2 coupled matrix evolution equa- 
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tions (3.85): 



r 



2 a f^"'^ + b r^-'^ + c { r^''^ [qq — App] + 2 r''''^ ^ [qP ~ P' 



[2pp-qq/2]+r^^+'-> [pq-qp]}/d, j = 1, 2, mod (2) 



^1) 



d 



2r ■ r 



q^ - 4p^ 



5 = 2 



1,(1) u;i) , u;2) 1,(2) 



^1) 



^2) 



^1) 



^2) 



(5.11a) 

(5.11b) 
(5.11c) 



This model is presented as a simple example of the kind of results that can be obtained 
by using the quadruplication (in fact, in this case, only duplication) technique described 
in Section |3.5| ; clearly many more models (which, however, become more and more com- 
plicated) can be obtained by iterated uses of these techniques, which can be moreover 
combined with those used above (complexification, and association with other solvable 
models, in particular to manufacture translation-invariant models). The explicit solution 
of this model (from (3.2), (3.84) and (4.3a)) is a straightforward, if tedious, task, that we 
leave as an exercise for the diligent reader. 

Next we consider the scalar /vector solvable one-body problem which is obtained by 
applying the parameterization (4.1) to the solvable matrix evolution equation (3.44). It is 
characterized by the following equations of motion: 



/9 = a + /3 p + 7 [p + c {p^ -r'^)] - c 



3 p p — 3 [f-f] + cp ip^ — 3r^) 



/3f + 7 



r + 2c pr 



3pf+3pf— rf\r + cr {3 p^ — r^) 



(5.12a) 
(5.12b) 



Note that in this case the 3-vector equation (5.12b) is coupled to the scalar equation 
(5.12a). One can of course complexify this model by setting, say, 

a = a + za, (3 = b + ib, 7 = c + ic, c = C + iC, (5.13) 
getting thereby a solvable 2-body problem. And one can perform an additional doubling 



via the technique of association., see Section |3^, getting thereby a solvable translation- 
invariant 4-body problem. 

Another solvable one-body problem is obtained by applying the parameterization (4.1) 
to the solvable matrix evolution equation (3.56). We do not report the corresponding 
equations, that the very diligent reader will easily get (using (D.8)). The remarks given 
above, after (5.12b), are then again applicable. 

Next, we report the equations of motion of the scalar/vector solvable one-body problem 
that is obtained by applying the parameterization (4.1) to the solvable matrix evolution 
equation (3.64): 



c+(a — a) p -\- aa p -\- ba (p^ 



2b [pp 



{a — a) r + aar + 2bapf — 2b {pf+pf \ —2bff\f 



(5.14a) 
(5.14b) 
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Note the similarity of these equations of motion to (5.12). The remarks made after that 
equation are as well applicable now. 

Application of the parameterization (4.1) to the integrable matrix evolution equation 
(3.50) (with the position Q = jl + iC ■ a) yields the following scalar/vector integrable 
one-body problem : 



2 c' 



2 c' 



— r r 



p {p'-^r')+^p-(c-r)' 

It is easily seen that these equations are yielded by the Hamiltonian 
ii-(p,7r;r,p) = (/-7r2)/2 + c2 

-27(rV)-4/0 (C.f) +r^-6r2/92 + /l /2. 



(5.15a) 
(5.15b) 



(5.15c) 



Of course only in the case with C = is rotation-invariance preserved. Then one can 
moreover set p = , so that (5.15a) is trivially satisfied, and (5.15b) reads 



af — br' f , 



(5.16a) 



where, for notational simplicity, we set c' = b/2 , 2c^7 = a . This is of course the one-body 
problem yielded by the standard Hamiltonian 



H{p,r} =p^/2 + V{r) 



with 



V{r) = -arV2 + 6r^/4 , 



(5.16b) 



(5.16c) 



whose intcgrability is well-known (the standard one-body Hamiltonian (5.16b) is of course 
integrable for any spherically symmetrical potential energy V {r)). Clearly all solutions 
of this model are confined if 6 > . An analysis of these motions is beyond the scope of 
the present paper, although we plan to take a look in the future at the interplay of the 
technique of solution entailed by the finding reported above, with the standard method to 
solve one-body spherically symmetrical problems in three-dimensional space. 

We leave as an exercise for the diligent reader the derivation of the equations of motion 
entailed by the application of the parameterization (4.3) to the integrable matrix evolution 
equation (3.62). 

Next we report the linearizable one-body problem that is obtained by applying the 
parameterization (4.3) to the integrable matrix equation (A. 8b) (using (D.IO)): 



r = 2ar -\- br-\- k 



2 r ( r • r 



/r' + 



r Ar . 



ir 



(5.17) 



Here / (r) is an arbitrary function (in fact, only its odd part contributes). If we restrict 
attention to the simplest case in which / (r) is linear in r, / (r) = Cr/2, and we moreover 
set A; = 0, we get: 



2ar + bf+ C r A r . 



(5.18) 
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Note that the above linearizable one-body problem can be directly obtained by applying the 
parameterization (4.3) to the matrix evolution equation (3.35) (with the trivial notational 
change c ^ C). A discussion of this problem and its solution can be found in Appendix 
B. 

Complexification of equation (5.18) via the positions 

f=f(^)+ifi'^\a = a + ia, b = f5 + ij3, C = c + ic, (5.19) 
yields the integrable 2-body problem characterized by the equations of motion 

r =2lar —ar j+pn'—p'n'+cl'n'Ar —'n'Ar 1 

- c (r-<i) A f^^^ + f=<2) A r ^^) , (5.20a) 



r^'^ = 2 (a r^'^ +af^'^)+P f<'^ + P r<'^ + c (r-<^) A r^'^ + r<'^ A f^'^) 

+ c (r-<i) A r^'^ - r-<2) A r^'^) . (5.20b) 

For an analysis of some aspects of the solution of this problem we refer to Appendix B, 
where in particular it is shown that, at least in the 2 cases characterized by the restrictions 
a = $ = 0, a = 3lo/2, (3 = 2 lo"^ or a = /3 = f3 = 0, a = u;/2 , with u> an arbitrary (real, 
nonvanishing) constant, this model is solvable and all its solutions are completely periodic 
with period T = 27r/uj. 

Association of (5.20) with an appropriate, trivially solvable, model of type (5.6) yields 
the linearizable translation-invariant 4-body problem characterized by the equations of 
motion 

r = ^ ^ s —7s ±[ 2 lar — ar \ -\- p-n ' — p-n > 

+ c (r-<i) A r^'^ - r-<2) A r^'^) - c {t^"-^ A r^'^ + r-<2) A r^'^) ] } /2 , (5.21a) 
r = \1 s +75 ib[2(ar + ar ] + pn ' + pn ' 



+ c (r-<i) A r^'^ + r-<2) A r^'^) + c {f^^ A r^'^ - 7^^^ A r^'^) ] } /2 , (5.21b) 

^1) = ^l>+) _ ^1 ^2) ^ ^2,+) _ ^2,-) ^ 

S-<1) = ^1,+) + ^1,-), s-<2) = ^2,+) ^ ^2,-) _ (5 21c) 

This model is solvable and all its solutions are completely periodic with period T = 27r/a; 
if there holds either one of the 2 sets of restrictions on the 4 coupling constants a, a,/3,/3 
reported above (after (5.20b)), and in addition there hold the 2 constraints 7 = 0, 7 = mu;, 
with m an arbitrary integer {m ^ 0; for m = 0, namely 7 = 7 = 0, (5.21) reduces to (5.20)). 
Last but not least, let us display a class of one-body problems, 

f{t) = (p{r)r{t)A'r{t) , (5.22) 

where (p (r) is an arbitrary function. The above equation can be obtained via the parame- 
terization (4.3) from the integrable matrix equation (3.36) (with a = 6 = 0, see Appendix 
A for details). The simplest choice (p{r) = k yields 

?{t) = k f{t) A fit) (5.23) 
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which is a subcase of (5.18) and has a remarkably neat physical interpretation (in spite 
of its parity- violating character): it describes the motion in ordinary space of a particle 
acted upon by a force proportional to its angular momentum (of course the arbitrary 
coupling constant k could be rescaled away). This is a solvable equation: its explicit 
solution in terms of parabolic cylinder functions is given in Appendix A (see (A. 19), (A. 28) 
and (A.36)). 

A solvable translation-invariant 2-body problem can be obtained from (5.23) via the 
positions (5.5b,c), with s{t) satisfying the solvable equation of motion (5.6b). Its equations 
of motion read {k = C): 

= {« (r^+) + r^-)) ± C [{f^+^ - f^-^) A (r^"'^ - r^"^)] } /2. (5.24) 

Next we consider the choice (in (5.11)) 

ip{r) = k/r'^, (5.25a) 
yielding 

m = . (5.25b) 

Also in this case we have an explicit solution in terms of hypergeometric functions (see 
Appendix A, formulas (A.19), (A.28) and (A.38)). 

A "physical" case obtains from (5.22) via the choice: 

ip{r) = k/r^ , (5.26a) 

yielding 

f (t) = UmAm , (5.26b) 

namely the "Newton / Lorentz" equation of motion of an electrical charge in the magnetic 
(radial) field of a magnetic monopole (or of a magnetic monopole in the electrical field 
of an electrical charge). This equation is partially investigated in Appendix A, a more 
complete investigation is postponed to a subsequent paper. 

We end this subsection with two linearizable three-body problems that are obtained from 
the "Nahm equations" (3.67). The first one is obtained by applying the parameterization 
(4.3) to the linearizable matrix equation (3.69) and it reads: 

•A 2 , 

Tn = |/^n+lMn+2r„+2 A r„+i 

3 1 
"I" ^ ^ [ '^(^n,mTm 
m=l 

+ i^n+l,ml^n+2rn+2 " an+2,mA'n+l^n+l + ^ an+l,man+2,fent^ A rv„]} . (5.27) 
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In the above equations as well as in the following ones up to the end of this Subsection 



5.1 , all indices are defined mod(3). 

The second linearizable three-body problem is obtained from (3.67) via duplication (see 
(3.80)) and by again applying the parameterization (4.3). The relevant formulas read: 



( An 






An 


i ^ 

— Uji ■ (7n 




C 




^ ^n ' ^n 





An = -Un ■ <7n , (5.28b) 

c 

= ifn ■ an , (5.28c) 

Their insertion in (3.67) yields, after some labor, the following Newtonian equations of 
motion: 

rn = [an + (? i/i,^! + ^n+2)] " [&n,n+l + (? (^n • ^n+l)] 

- [&n,n+2 + ^ {rn ' r„+2)] f„+2 , (5.29a) 



where: 



{ Un+i + + (7n,n - 7n+l,n+l " 7n+2,n+2) [(r*n+l ' Un+l) + (^1+2 • Un+2)] 

- 7n,n+l [{rn ' "^n+l) + (f„+i • Un)] - 7„,„+2 [(rk • Un+2) + (r„+2 • Un)] } / , 

(5.29b) 

bn.m = { Un-Um + ^ {lm,m - 7m+l,m+l " 7m+2,m+2) [(^n • Urn) + (?V« " ««)] 
7m, n 

[(f„+i • Un+i) + (f„+2 • ?i„,+2)] 
+ 7m ) + (r2n-m • Un)] } /A^ , (5.29c) 

7n,m = rn-fm , (5.29d) 

Un = (7n,n " 7n+l,n+l " 7n+2,n+2) fn + 7n+l,nni+l + 7n+2,n^\^+2 , (5.29e) 

A = fi • fa A fg . (5.29f) 

Note that an apparent paradox arises, namely, setting c = in the above formulas, the 
"free motion" equation f = is not recovered (in contrast to what seems suggested by 
(3.67)): this is due to the definition (5.28b), where the coupling constant c appears in the 
denominator. 

5.2 Many-body problems 

In this subsection we consider several solvable and/or integrable and/or linearizable many- 



body problems, obtained by applying the (first) multiplication technique of Section 3.5 



and Appendix C to some of the models of the preceding Section 5^, or some of the 
parameterizations of Section ^ to some of the matrix ODEs of Section |^. 

The first model we report is obtained by applying the first multiplication technique of 



Section 3.5, or Appendix C, to the solvable equation of motion (5.18). 



The relevant equations of motion then read 

N N 

^2 dn— ni rni ~\~ ^n— ni '"ni^ ~l" ^ ^ ^Cn— m— n2 '"ni A ?'n2^ • (5.30) 



ni=l ni,n2=l 
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N 

These equations of motion are obtained from (5.18) via the positions r = Vn/Tn, a = 

n=l 

N 

Yl Vn flnaiid SO on (see Section 3.5, or Appendix C), which entail that the indices n, ni, n2, 

n=l 

see (5.30), run from 1 to and are defined mod(A^). For instance for = 3 these 
equations of motion, when displayed in longhand, read 

fi = 2 ( ai fs + 02 r2 + as fi ) + bi + 62 ^2 + ^3 n 



Cl r*3 A fs + C2 n A fi + C3 f2 A f2 



+ 2 



+ Cl ri A r2 + r2 A ri + C2 ?"2 A ra + r2 A ra + C3 [ri Ar^ + r^ Ari 



(5.30a) 



f 2 = 2 ( ai fi + 02 rs + as f2 ) + bi n + 62 ^3 + ^3 r2 



Cl f2 A f2 + C2 rs A fs + C3 fi A fi 



+ 2 



+ci ri A r3 + rs A ri + C2 ri A r2 + ri A r2 + C3 r2 A rs + rs A r2 



(5.30b) 



fs = 2 ( oi f2 + 02 fi + as fs ) + 61 f2 + 62 fi + ^3 r3 



Cl ri A ri + C2 ^2 A r2 + Cs rs A rs 



+ 2 



+ Cl A rs + r2 A rs j + C2 A rs + rs A ri j + cs A r2 + r2 A ri 

(5.30c) 

Let us however recall that the equations obtained via this technique can generally be 
uncoupled via a linear transformation (see (3.92)). 

Association of the model (5.30) with the trivially solvable model characterized by the 
equations of motion (see (5.6)) 



C^n ^r. 



(5.31a) 



yields, via the positions. 



1"), rl±) = (s„±r-;)/2. 



(5.31b) 



the following equations of motion characterizing a linearizable translation-invariant many- 
body problem: 



■4±) 
r 



N N 

{ Ckfi Sn i ^ ^ ^2 ajj_m r^^ + bn—m ^ni^ ~l~ ^ ^ (^Cn—ni~n2 ^ni A ?'n2^ } /2 • 
ni=l ni,n2=l 

(5.31c) 



In these equations of motion for the 2 A^ three- vectors rl^^ , the three- vectors s„ and 
f„ are defined by (5.31b). Of course all indices are again defined mod(A^). 
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Likewise, let us mention the solvable Hamiltonian many-body problem characterized by 
the equations of motion 

N N 

Qn ~ ^ ^ (o^n— ni Pni "I" Tn— ni Qni) + ^ ^ [Cn— ni— n2— ns iPni ^ Qn2^ ^ Qnz\ t 
ni=l n\,n2,nz=l 



N N 

Pn ~ ~ ^ ^ iPn—ni Qni "I" Tn— ni Pm ) + ^ ^ [C; 
ni=l ni,n2,n3=l 



(5.32a) 

n— m— n2— ns (Pm ■''^ 9712) ■''^ Pns] j 

(5.32b) 



whose derivation is relegated to the end of Appendix F. But we exhibit here in longhand 
the form these equations of motion take for AT = 2: 

=aip2 + a2 pi + 71 ^2 + 72 

+ ci [{pi A ^i) A ^2 + {Pi A ^2) A gi + {p2 A qi) Aqi + {p2 A 5*2) A 5*2] 

+ C2 [(pi A qi) A gi + (pi A ^2) A g2 + (P2 A gi) A ^2 + (#2 A ^2) A gi] , (5.33a) 

92 = aiPi + Q;2P2 + 7i 91 + 72 02 

+ ci [{pi A gi) A gi + (pi A g2) A g2 + {p2 A gi) A g2 + (p2 A g2) A gi] 

+ C2 [(pi A gi) A g2 + (pi A g2) A gi + (p2 A gi) A gi + (p2 A g2) A g2] . (5.33b) 

Next, we report the scalar/vector solvable many-body problem characterized by the 
equations of motion 

N 

Pn ^" ^ ^ {Pn—ni Pm ~^ In—ni Pm } 

ni=l 
N 

3 ^ ^ -^Cj^—^j—jij Pm Pn2 {j^m ' ^ 

rai,ra2=l 

N 

"I" ^ ^ {Cn— m— 712— ns 7ni [Pn2 Pns (''n2 ' ^ns)]} 
ni,n2,n3=l 
AT 



')]} 



E {< 



-n— ni— n2 



-ns—Ui Cni Pn2 [Pns Pn4 3 (j'ns ' '"n4)]} j (5.34a) 



ni,n2,n3,n4=l 

N 

"n ~ ^ ] "l/^n— ni ''ni 7n— m ^ni 



ni=l 

AT 

^ ^ "Icfi—yij^—fij Pm ^n2 "I" 3 Pni '"n2 {^m ^ '"n2^ j ^ 
ni, 712=1 

N 

"I" 2 ^ ^ {Cn— ni— n2— ns 7ni Pn2 '"nsi 

ni,n2,n3=l 
AT 

~ ^ ^ {Cn— ni— n2— ns— n4 Cni ''n2 Pns Pn4 ~ (^ns ' ^n4)]} > 
ni,n2,ns,n4=l 



(5.34b) 
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which are obtained by applying to (5.12) the (first) multiplication technique described 
at the end of Section |3]^ and in Appendix C. Hence these equations of motion can be 
uncoupled by using a linear transformation of type (3.92). 

Next, we report the equations of motion of the solvable many-body problem which is 
obtained by applying the parameterization (4.3) to the special case of the solvable matrix 
evolution equation (3.63) with D. = D^ = X (see Appendix H): 



[rn ■ Tn] fn + ttn 7n A fr. 



fn = anfn- a-n In ^'n A f „ + { - 7„ f „ A 
N 

+ lnrn/\ [7n (^n ' ^m) Tn + a\ Tm + 7„ fn A rm]}}} 



m=l 



(5.35) 



Here the (3 + A^) constants a„, a^; 7™, bnm are arbitrary. An analysis of the motions 
entailed by this model is straightforward (see Appendix H), but it exceeds the scope of 
this paper. Note that even the one-body case (A^ = 1; or, equivalently, bnm = ^nmbn) is 
nontrivial (namely, its equation of motion is nontrivial; of course in some sense all solvable 
models are trivial!). 

Next, we report the equations of motion characterizing the scalar /vector solvable many- 
body problem which is obtained by applying the parameterization (4.1) to the solvable 
matrix evolution equation (3.65): 

Pn = {an - CLn+l) Cn + {an - a„+i) (a„ - Qn) Pn + Cn {bn Pn " ^n+l Pn+l) 

- {an - an) bn+1 {pn Pn+l - Tn ■ f„+i) - 3 6„ [pn Pn - Tn ■ 
+ {dn-2an + ffln+l) [Pn + bn {Pn " 

+ bn+l {pn Pn+l " ' ^n+l) + K bn+1 [Pn+1 {Pn - ^1) - 2 /9„ f„ • fn+l] 

- bl Pn {pi - 3r^) , (5.36a) 

rn = {an - an+l) {an - an) fn + Cn {bn r„ - fn+l) 
bn+1 {an an ) {pn fn+l + Pn+l fn + fn^ fn+l) 

bn ^3 Pn fn + ^Pnfn-fnf\ fr^j + (a„ - 2 a„ + a„+i) {fn + 2bn Pn fn , 
+ bn+1 ( Pn+l fn + Pn fn + rn A rn+i 



+ bn bn+1 [2 Pn+l Pn fn + (Pn " ^1) fn+l +2pnfn^ fn+l] " bl (3 pi - r^) f„ . 

(5.36b) 



Next, we report the equations of motion that are obtained by applying the first tech- 
nique of multiplication of Section 3.5 and Appendix C to the integrable one-body problem 

N 

(5.15) (with C = 0, to restore rotation-invariance, and via the positions 2c^ = ^ ??nbn, 

n=l 
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N N N 

2 7 = ^ a„ as well of course as p = ^ 'f]nPm^= ^ T^n^n)- They read: 

n=l n=l n=l 

N N 

Pn — ^ ^ {O'n—ni Pn\) + ^ ^ {^n— ni— n2— na Pni [Pn2 Pna ~ 3 (r^j " ^na)]} j 
ni=l n\,n2,nz=l 

(5.37a) 

Af TV 

'"n — ^ ^ ('^n— ni ^ni ) ~l~ ^ ^ {^n— ni— n2— na ''ni [3/5^2 Pna ('"^^2 ' '"'^s)]} • 
n\=l n\,n2,nz=l 

(5.37b) 

In these equations all indices are defined mod(A^). This scalar/vector many-body prob- 
lem is then of course an integrable one, but susceptible to being uncoupled via a linear 
transformation of type (3.92). 

We plan to discuss in a separate paper the integrable class of (rotation-invariant!) 3- 
dimensional many-body problems with quadratic and quartic interactions obtainable by 
applying to the integrable matrix equation (3.50) other techniques of association and 
multiplication (see Section |3.5| ) as well as appropriate parameterizations of matrices in 
terms of 3- vectors (see Section ^) . 

Next we report the linearizable many-body problem whose equations of motion are ob- 
tained by applying the parameterization (4.3) to the linearizable matrix evolution equa- 
tions (3.29), and therefore read as follows: 

N 

f„ = ^ {2 a„„, (t) + bnn, (t) f„, + Cn, (t) 2 (^f„,, A 

ni=l 

- Cni {t) (f„ Af„J} 
N 

+ {2anm (t) {t) {rni A - Cm (t) C„2 (i) [{fn A r„i) A f^^]} . 

ni,n2=l 

(5.38) 

The connection with the notation of (3.29) is: C„ = 2c„. Note that, as in (3.29), we are 
allowing in this case all the N {2N + 1) quantities anm (t) , bnm {t) , Cn (t) to be time- 
dependent. 

We leave as an exercise for the diligent reader to derive the more general (and of 
course also linearizable) equations of motion that are obtained by applying to (3.29) the 
parameterization (4.1) (rather than (4.3)). 

Next we report the integrable many-body problem whose equations of motion are ob- 
tained by applying the parameterization (4.3) to (3.38), leaving again as an exercise for 
the diligent reader to derive the more general (and of course also integrable) equations of 
motion that are obtained by applying to (3.38) the parameterization (4.1) (rather than 
(4.3)). They read: 



2 r^ I r^i ■ r^i I I r^ ■ r^ 



I n 



+7 {^n+i - [2f„ (f„ • f„_i) - f„_i (f„ • f„)] /r^_i} . (5.39) 



356 



M. Bruschi and F. Calogero 



Note the "nearest-neighbour" character of these equations of motion (we do not repeat 
here, nor below in analogous cases, the comments on the need that these equations of 
motion be supplemented by "boundary conditions at the n-ends:" see Sections ^ and 3.3). 

Another scalar/vector integrable many-body problem is obtained by applying the pa- 
rameterization (4.1) to (3.39b) (note that in this case the parameterization (4.3) is not 
applicable; on the other hand the simple, quadratic, character of the nonlinearity, see 
below, makes the following equations of motion particularly suitable for the multiplica- 
tion techniques of Section 3^; but we leave this additional extension in the many-body 
direction as an exercise for the diligent reader). The equations of motion read 

Pn = c pn (pn+i - 2 /j„ + p„_i) - f„ • (f„,+i - 2 f „ + f„_i) , (5.40a) 
rn = c [fn {pn+1 - 2 /J„ + p„_i) + (f„+i - 2 f„ + r„_i) 

+ rnA {fn+l - fn-l)] ■ (5.40b) 

Note the translation-invariant character of these equations of motion (i. e., invariant 
under Pn ^ Pn + Po ^ ^ + fo with pq, ff) time-independent but otherwise arbitrary). 

Let us also report the generalization of these equations that are obtained by first re- 
placing formally in (5.40) p„ (t) , fn (t) with, say, pn (r) , fn {t) (also replacing, of course, 
derivatives with respect to t with derivatives with respect to r), then setting Pn{T) = 
exp (—at) Pn (t), fn (r) = exp {—at) fn (t), t = [exp (at) — 1] /a and then also performing 
a standard complexification (p„ = pn^ + i pn^ , f„ = f^^^ + i rl^^ , = 7 + ^7 , a = X-\-iuj). 
We thus obtain the more general scalar/vector integrable many-body problem characterized 
by the following equations of motion: 

+ 7^^ -7i„ + (A7-a;7) 5„ - (A7 + W7) Bn , (5.41a) 

+ 7A„ + 7i„ + (A7 + w7) Bn + {X-r-^l) Bn , (5.41b) 

+ 7C„-7^„ + (A7-(J7) ^„-(A7 + w7) 5„ , (5.41c) 

+ 7(7„ + 7d„ + (A7 + (J7) 5„ + (A7-w7) 5„ , (5.41d) 

A = 4(1) _ 4(2) A =4(1) ,4(2). R = r(1) _ r{2) r =r(1) I r{2) 
-^n — -^n ' -^n — -^n ' ' — J^n ' — J^n ' i 

(5.41e) 

aU) = A(i) (Si) _ o Si) _L Si) \ I 'S^) . ( di) _oSj) 1 di) ^ i - T 9 

(5.41f) 

Aii) = (o^^^ - 2 o^i) + o^^'^ ) + f^^"^^^ • (f^^'^ - 2f^i) + f^^^ ) 

^n — Pn I Pn+1 ^ Pn ^Pn-lJ^'n l'n+1 ^'n ^'n-ll ' 

j = l, 2, mod (2) , (5.41g) 

nii) = Si) (Si) _ 2 Si) + S^^ ^ + f<i) ■ ( f^^^ -2r^i) + f^^'^ ^ 1-12 

-°n — Pn \^Pn+l ^ Pn Pn-lj 'n y n+1 ^'n ^'n-lj^ J — ■•-) ^ ; 

(5.41h) 

-Dn — Pn \^Pn+l ^ Pn Pn~l^ ' n \^'n+l ^ 'n ^ ' n~l ) ' 

j = 1, 2, mod (2) , (5.41i) 



Pn 


-A 


Pn 




Pn 


-A 


Pn 




•41) 

'^n = 


-A 


ss 

^n 


^ 42) 


42) 
'^n = 


-A 


42) 
'^n 


41) 



'^n — ^n i '^n — '^n ^ ^n i — ^n ^n i ^r, 



n 5 



(5.41j) 
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(5.41k) 



c = f* 



A 



flV^^ _ 9 rt'-''' -I- n^-^^ -I- h^J 
Pn+1 Pn ^ Pn-1 ) ^ Pn 



0+1) 



(r-gi-r-J2i) , i = l,2, mod (2) , 

—'n \Pn+l Pn ^Pn-lJ^Pn y n+1 ^'n ^'n-lj 



n+l n— 1 



, J = 1,2 



— 'n \Pn+l ^ Pn ^ Pn-lJ ^ Pn y n+l ^ ' n ^'n-lj 

+ z r-<i+i) A f ri^l, - /Jl, ) , .7 = 1,2, mod (2) . 



n+l ' n-1 I ' 



(5.411) 
(5.41m) 

(5.41n) 



Note that, as implied by the derivation of these equations of motion, we can assert that 
all their solutions are completely periodic, with period T = 27r/a; , if A = 0, co ^ 0. 

Next, we report the scalar/vector linearizable many-body problem characterized by the 
equations of motion 



Pn = [pI Pn + 2pn (jn ' Tn) - Pn (r„ ■ r„) / [pi + r^) 

"I" ^ ^ {^nm Pm Pm Pn ~^ Pm (Xm ' ''n) Pm (^m ' 1"nj 

m 

m ■ fm)]/ {pm + ^m)} i 
rn = {rn Pn Pn + (^n " 4^ - 4 Pn+ {fn ' ^n) 



(5.42a) 



rn [p 



+ fn Pn Pn 



}/{pl + rl) 



Pm Pn + 



Pm Pm ~\~ I fm ' f; 



+ 



Pm ^m ^ Pm fm A Tn ~\~ Pn ^m ^ 



}/{pm + ri) 



(5.42b) 



which are obtained by applying the parameterization (4.1) to the linearizable matrix ODE 
(3.47). 

Next, we report the linearizable scalar-vector many-body problem that is obtained by 
applying the parameterization (4.1) to the linearizable matrix ODE (3.43d). Its equations 
of motion read 

Pn = apn + PnPn - Tn ' Tn + C ^Pn-lPn - fn-1 ' "Tn - PnPn+1 + ^n " ^n+l} , (5.43a) 
fn = fn "I" Pnfn "I" Pn^n "^n ^ fn 

n fn Afn+i} , (5.43b) 



with 



Pn = [pnPn + • / [pi + r^) , 



(5.43c) 
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rn = {pnTn - PnTn - A / (p^ + r^) , (5.43d) 

and 

Pn+l = PnPn+l " ' r„+i , (5.43e) 

= PnTn+l + Pn+l^n - A r„+i . (5.43f) 

We do not write in explicit form the equations obtainable from this one by complexification; 
their interest is related to the fact that, if in (5.17a) o = ibia;, a; > 0, its generic solution 
is presumably completely periodic, with period T = 2t:/u. 

The diligent reader may try and derive the more general models that are obtained 
by employing the parameterization (4.6), which is applicable to (3.43d) (see (4.6m)) and 
yields equations of motion only involving 3-vectors, but in a more complicated manner. 

Next, we report the equations of motion of the linearizable scalar-vector many-body 
problem that is obtained by applying the parameterization (4.1) to the linearizable matrix 
ODE (3.66): 



Pn 



(On+l - On - bn) Pn + PnPn - fn ' Tn + &n+l (^Pn+lPn - • fnj 

+ Cn+lPn - Cn {pnPn-l " " , (5.44a) 

Tn = (On+l - On - 6n) + PnTn + PvTu -Tnl^Tn 

+ hn+\ (pn+irn + Pn^n+i - f„+i A f„ 



+ c„+ir„ - Cn (^Pnrn-1 + Pn-\^n + A f„_i^ , (5.44b) 
where Pn-,fn are given again by (5.43c), (5.43d) and 

Pn-l = PnPn-l " " Tn-l , (5.44c) 

rn-1 = PnTn-l + ?n-l^n - A , (5.44d) 

'Pn = {pnPn + r„ • f„^ / {p^ + " ^n) , (5.44e) 

Tn = {pnTn " PnTn + A / {p^ + ' ^n) • (5.44f) 

The highly nonlinear character of these equations of motion should be noted.We leave as 
exercise for the diligent reader to write the (only 3- vector but more complicated) equations 
of motion that obtain by using the parameterization (4.6), instead of (4.1). 

Next, we display the scalar/vector linearizable many-body problem that is obtained by 
applying the parameterization (4.3) and (4.1) to the linearizable matrix ODEs (3.23), via 
the positions C/„ = ifn-^,F = if-a,G = iq-d,H = ih-d,V = rj-\-iiI-d,Y = 9-\-iy-a: 



+ E {fm-2kArm}+ 6^"^ fm 

m=l 



P"n^n = 4^k Ar'n + '^k Afn- 4:k A{k A f„) 
N 
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+ a^r^) {/- 2 A /} + /+ 4"'^ - 2 A 5} 

+ 6^"^) g + a^,"'*) {K-2^A^} + 6^"'^^ /I , (5.45a) 
/x(^)/=4^a/+2^A/-4^a(^A/) + ^ [a(^«) ^fm-2k ^fm} +b\l''^ f, 



m=l 



+ a(-^^) |/- 2 A /} + b^ff^ f + a^f3) 2 A 3} 

+ h^^3) ^ + ^^h-2k^^+ h^f^^ h , (5.45b) 

//(^^) ^=4feA^+2^A5-4feA(feA5) + i; [a^^") {f^ - 2 A f^} + b^X^rm 

m=l 

+ a^^-'^) | / - 2 A /} + b^s-f^ / + a(»s) |^ - 2 A 5} 
H^^^ h = 4k Ah + 2k Ah - 4k A (kAh^ + ^ [a^/j") {fvn - 2 A f^} + fe^^") 

m=l 

+ a(''-'') |/-2fc A/} + 6('^-^)/ + a('*») |^-2fcA5} 
uS'"'>'i= 2kAi) + kAv + k'^v 



(5.45d) 



(vu) 2 
m 



-[^-kAv] a(^-) {r'^-2fcAf^} + 5^^ 

* 2 ^ 

+ a^''-^) {/ - 2 A /} + f + a^''^) {^-2^A^} + / 



+ a(^^) {/ - 2fc A Z}^ + b^'^-^) /2 + a(^^) - 2^ A5}' + b^^^) g 
+ 0^"^^ {^fl-2kAhy + b^''''^h'^} , 



,ivu) 



{n„-2^Ar^} +6(^")r; 



(5.45e) 



^Wy'= -2k Ay - k Ay + k"^ y 

N 



-[^+kAy}{Y^ \a^Z'^ {rm-2kAfm} + 6^"^ r 



m=l 
^1 2 



+ a(^^) { / - 2 A /} + 6(2//) /2 + ^{ra) |^_2fcA^} + 5^ 

2 ^ r 2 



m=l 
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2 2 

+ a^yf^ {/ - 2 A /} + 6^^^) f + a^y^^ [p-2k^g^ 

+ b^y<^^ 5^ + a^y^^ {K-2fcAK}V b^y^^ } , (5.45f ) 

-2 ?) — ^fc • i/j — ?7 — ^fc • i/j + (A;^ — 7^) -q 

N r 2 

-^{E {r''m-2^Ar^} +6(:")r^ 

* 2 ^ 

+ a^''-^) {/ - 2 A /} + /2 + a^'^y^ {^-2^A^} + t^"^) 5^ 

2 ^ r 2 

+ a('''^) {K-2fcAK} ^ a(;^«) {^^-2^Arv„} + ft^^^V, 



2 2 



N p 2 



(5.45g) 



* 2 2 

+ 0(2//) I / - 2 A /} + b^yf^ f + a(f ||'-2^A5'} + Uyy^ 



m= 1 

2 ^ 

+ 0(2//) { / - 2 ^ A /} + 6^^^) f + a^yy^ ^g-2k^g^ + fe^^/ff) ^ 
+ a^y^^ |K - 2fc A /i 



2 



2fc Arv 



(5.45h) 



where 



k = f + t]9 + (^f^^9^f + f^^y^ 1 = 9 ■ V + h ■ y . (5.45i) 

These equations of motion determine the evolution of the iV + 5 three- vectors Vn (i), 
/(t), ^(i), andofthe2scalarsr/(i), e{t). They feature (2 iV^ +20 iV+42) 

arbitrary "coupHng constants," hence they inchide many special cases, corresponding to 
appropriate choices of these coupling constants, many of which could be set to zero to 
obtain simpler systems; but the exploration of these cases is beyond the scope of this 
paper. 

The diligent reader may write the more general equations of motion involving + 5 
scalar quantities, as well as three- vectors, which obtain by using the parameterization 
(4.1) for all the + 5 matrices that evolve according to (3.66), rather than for only those 
2, y and Y_ , for which a parameterization (4.1), rather than (4.3), is mandated by the 
very structure of (3.23). 

Finally we display the equations of motion of the scalar /vector solvable many-body 
problem that is obtained by applying the parameterization (4.1) to the matrix ODEs 
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(3.96): 



N N 
+ 

mi=l mi,m2=l 



Pnmi Pm2in ( ^nmi ' ^m2m ) ^ Pnmi P'm2in 2 ^ Tnmi ' ^m2m , 
N 

~l~ ^ ^ '^mim,2 ^msm [Pnmi Pm2m3 i^nmi ' ^7712^3)] 
mi,m2,m3=l 

AT 

^ ^ dmim2 drnznii [Pnmi Pm2mz Pnum Pnmi {^m2m3 ' Tm4,m) 

mi,m2,r?i3,r?i4=l 

Pm2m3 (^nmi ' '^rnum) Pm4m iXnini ' '^m2mi^ ~^ (^nmi ^ ^m2m3) ' ^m4m ] j 

(5.46a) 

AT 



mi=l 

AT 



Tnm — / ^ Y^mim ' nrrii ~r t-mim ' nmi 
mi=l 

AT 

m2 fnm,\ ^ fm2m ~\~ 2 fnmi ^ ^m2 



mi,m2=l 
AT 

^ ^ dmim2 ^msm fnm\ ^ ^m2mz 
mi,m2,m3=l 
N 

^ ^ C^m3m4 f^mim2 

mi ,m2,m3, 7714=1 

[ ^nmi {Pm2mi Pm^ra (fm2mi ' ^7714771)} ~l~ ^77127713 {Pnmi Pm^m iXrimi ' '"7714777)} 
+ fm4m {Pnmi Pm2m3 {l"nmi ' '"77127713 )} ] • (5.46b) 

These equations of motion describe the evolution of the A^^ scalars pnm {t) and of the 
N'^ three-vectors fnm (t); they feature 4 A^^ arbitrary "coupling constants." 



6 Outlook 

This paper may have appeared extremely long to the researcher, if any, who had the 
stamina to read through it; yet it is clear that here we merely introduced our topic. Much 
remains to be done, in the way of more systematic exploration of solvable and/or integrable 
and/or linearizable many-body problems in 3-dimensional space, as well as analysis of their 
detailed behavior and identification of interesting applications. Indeed, as emphasized in 
Section ||, the purpose and scope of this paper has been to introduce various techniques 
for the identification of such models, rather than to explore completely their potentialities. 
We are moreover aware of other developments [10,11] which might well provide additional 
starting points for applying the techniques described herein but which have not been 
discussed above to avoid our paper becoming excessively long. 

Let us end this section by indicating three additional developments beyond those we 
just mentioned, which are also suggested as a follow-up to our findings. One is the de- 



362 



M. Bruschi and F. Calogero 



tailed study of all cases in which the many-body Newtonian equations of motion exhibited 
above, as well as others that may be obtained by analogous techniques, can be fitted into a 
Hamiltonian framework; their investigation is important in itself, given the special interest 
of Hamiltonian systems, and even more so inasmuch as it opens the way to quantization, 
with the expectation that many-body problems treatable in a classical context be also 
amenable to exact analysis in the quantum case. Secondly, a vast vista of generalizations 
is implied by the possibility to consider discretized problems (in space and/or time). And 
thirdly, there is the possibility to manufacture many-body problems which feature New- 
tonian equations of motion with velocity-independent forces, and are amenable to exact 
treatments (solvable, integrable, linearizable) only for a subset of initial conditions (typi- 
cally, with arbitrary initial positions, but with initial velocities determined by the initial 
positions): these are the systems characterized by first-order equations of motion, in which 
the second-order ( "Newtonian" ) equations of motion are obtained by time-differentiating 
the first-order equations and then, in the second-order equations, using the first-order 
equations equations to eliminate the (once-differentiated) terms. 



Appendix A. On the matrix ODE U_ = 



LL,fiU) 



In this appendix we mainly focus on (2 x 2)-matrices, and in this context we treat (indeed, 
in some cases, solve) the matrix ODE displayed in the title. But first we tersely discuss 
the theoretical background to the linearizability of this matrix ODE, and indeed more 



general versions of it, without any restriction on its rank (see Section 3.2 for an analogous, 
albeit more general, treatment). 

Let the (time-dependent) matrices M(t)and U_{t) be related by a (time-dependent) 
similarity transformation: 

u{t) = W{t)U{t)[W{t)]-\ U{t) = [Wit)]-^u{t)Wit), (A.l) 

with the time-evolution of the (invertible!) matrix W_ (t) characterized by the matrix ODE 

W{t) = W{t) l{U{t)) , (A.2a) 
W{t) = l{u{t)) W{t) . (A.2b) 

The choice of the function / (U) remains our privilege; but note that, as entailed by 
the consistency of (A. 2a) with (A. 2b) via (A.l), we assume f {U_) to be a scalar/matrix 
function (i.e., a function depending on scalar quantities and on the matrix U_ , but on no 
other matrix, so that / (u) = Wl{U) W'^, [{U) = W'^ [{u) W, see (A.l) ). 

The matrix evolution equation (A. 2) defines the matrix W_{t) up to the assignment of 
an initial condition, 



W{0) = Wo, (A.3) 

whose choice remains our privilege, except for the requirement that the matrix W_q be 

invertible, see (A.l). Generally it will be convenient to make the very simple choice 
Wq = 1 (see below). 
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It is now easily seen that (A. 2a) entails 
u = WUWr\ U=W'^uW, 



u = w[u+[l{u),u]]w~\ u+[l{u),u 



(A.4) 
(A.5) 



Here and throughout [^A, fi] denotes the commutator of the two matrices A and 5, 
B\ = AR - BA. 

Assume now the matrix u (t)to evolve in time according to a second-order linearizable 
ODE, say 



u (n, u) 



(A.6) 



Here the choice of the function u(u, n) remains again our privilege, subject to the re- 
striction that (A.6) be linearizable (or perhaps even solvable] indeed the case on which 
we will mainly focus below corresponds to the simply choice u{u,u) = 0); but we must 
again assume that u is a scalar /matrix function of its 2 matrix arguments (namely that 
it depends on scalars and on the 2 matrices uand u, but on no other matrices, so that 
u (u, ii) W_ = u (jZ^Llji see (A.l) and (A.4)). It is then clear that the corresponding 
evolution of L[(t), characterized by the matrix evolution equation (implied by (A.6) via 
(A.5), (A.l) and (A.4)) 



U = u U,U] + 



uJiu) 



(A.7) 



is linearizable. Indeed its solution can be obtained via the following steps: (i) given 
the initial data U_{0) , C^(0), evaluate, via (A.l) and (A.4) with (A. 3), the initial data 
(0) , u (0) (note that the convenient choice W_q = 1 entails simply n (0) = U_ (0) , u{0) = 
U_{0) ); (ii) evaluate u{t) from the (by assumption linearizable) evolution equation (A.6); 
(iii) evaluate W_ (t) by integrating the linear nonautonomous matrix ODE (A. 2b) (with 
(A. 3)); (iv) evaluate lJ_{t) from the second of the (A.l). 

For instance by assuming that u{t) satisfies the solvable matrix ODE (3.1), namely by 
setting in (A.6) (up to trivial notational changes): 



2 ail + bu + kiiu 



-1 



(A.7) becomes 



U = 2aU + bU + ktJU~^U + 



u,fm 



(A8a) 



(A.8b) 



and more specifically (to make explicit the connection with the treatment of Section |3.2| 
by moreover assuming A; = , f {U_) = cU so that u(t) satisfies the solvable (linear! 
matrix ODE (3.37a), namely: 



u = 2au + bu , 
(A.7) becomes (3.36), namely: 



(A.9a) 



U = 2aU + bU + c 



U, U 



(A.9b) 
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The solution of equation (A. 9a) (with a and h two arbitrary constants) reads 

u (t) = exp (a t) [u (0) cosh {At)+u (0) A'^sinh (A t)] , (A.lOa) 

A= (a2 + 6)^/^ . (A. 10b) 

Hence the corresponding evolution equation for W_, see (A. 2b), reads 

W = e^p{at) [^cosh(Ai) + 5A-^sinh(At)] W, (A.ll) 

with A = cu{0), S = cu{0) two constant matrices (if also the scalar c is a constant). 

In Appendix B an additional analysis is given of the linearizable matrix evolution 
equation (A. 9b), and in particular it is shown that, if a = (3/2) iu, 6 = 2a;^ , or if 
a = iuj/2, 6 = 0, with lo an arbitrary (nonvanishing) real constant, (A. 9b) is in fact 
solvable and all its solutions are completely periodic, with period T = Itt/lo. 

The exploration of additional versions of (A. 7), corresponding to other choices of the 

(scalar /matrix) function u(^U_,Uj and f{U_), is left as an amusing task for the diligent 
reader. 

Hereafter in this appendix we focus on the simple choice 

f{u) = ^h{t)u, (A.12a) 

and u = 0, entailing (see (A. 6)) 

u = , (A.12b) 

hence 

u{t) = u{0) + tu{0) . (A.12c) 
Then (A.7) reads: 

U = ^h{t) U,U . (A.12d) 

Hereafter our attention is restricted to (2 x 2) matrices. Indeed our main aim now is 
to solve the 3- vector Newtonian equation of motion (corresponding to (A.12d) via (4.3)) 

■F{t) = h (r^) f{t) A f{t) , (A.13) 

where the function h (r^) is defined below, see (A. 15). 

This evolution equation clearly entails r-f=f-f=0 , from which we easily get 

f (t) = b , (A.14a) 
r^{t) = \f{t)f = x{t) = a + 2pt + bH^ , (A.14b) 
rit)-rit)=P + bH = ^icit) , (A.Mc) 
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with the 3 constants a, /3 and b defined in terms of the initial data as follows: 



a 



[f(0)r , /3 = f (0) • f (0) , b= f(0) 



(A.14d) 



{Beware: the quantity b introduced here, see (A. 14a) and (A.14d), has no relation to 
the "coupling constant" b in (A. 8) or (A.9)). Taking into account (A. 14b), it is clear that 
(A. 13) is obtained from (A.12d) via (4.3), with 

h{t) = h{x{t)) = h{r^) . (A.15) 

If we now time-differentiate (A. 13), taking into account (A. 14b) and (A.15), we get 



h/h) f+ 



(A.16) 



We thereby see, via (A. 14), that each component of the 3- vector r satisfies the following 
third-order linear nonautonomous ODE (written for the dependent- variable F (t)): 



F{t) = h{t)/h{t) F{t) + h^{t) 



1 



x{t)F{t)-x{t)F{t) 



(A.17) 



Let us define 3 independent solutions, Fq,Fi,F2., of the above ODE, via the following 
initial conditions: 

Fo(0) = l, Fo(0) = 0, Fo(0) = 0, (A.18a) 

Fi(0) = 0, Fi(0) = l, Fi(0) = 0, (A.lSb) 

F2 (0) = 0, F2 (0) = 0, F2 (0) = h{x (0)) = h{a) = h (0) . (A.18c) 

It is then easily seen that the solution of (A. 13) is given by the formula 

r{t) = f (0) Fo (t) + f='(0) Fi {t) + f (0) A f='(0) F2 (t) . (A.19) 

The task to solve (A.17) directly is daunting: nevertheless formula (A.19) provides an 
useful representation of the solution r{t) of (A. 13) in terms of the functions Fq, Fi , F2. In 
order to find them, instead of solving directly (A.17) it is convenient to use the technique 
outlined at the beginning of this appendix (and, more generally, in Section 2^) . Let us 
summarize the procedure: the solution of the 3- vector equation (A. 13) is obtained from 
the solution of the matrix equation (A.12d) (using the representation (4.3)); the matrix 
solution of (A.12d) is given by (see (A.l)) 



u{t) = mt)]-'u{t) w{t) 

with (see (A. 2a), (Al2a)) 

W{t) = \h{t)u{t)W{t) , 
W{0) = 1, 
and (see (A12.c) and (4.3), (A. 20a)): 



(A.20a) 

(A.20b) 
(A.20c) 



u{t) = i f(0) + f(0) t 



a 



(A.20d) 
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Hence we must now solve the (2 x 2)-matrix equation (A. 20b) with 
a + bt c 



uit) 



(a + bt) 



(A.21a) 



which corresponds to (A.20d) via the convenient choice we make for the cartesian reference 
system, so that its z-axis point in the direction of the 3-vector r= (a + 2/?t + 6^t^) = 

[x(i)]^''^, and its xz-plane contain the 3-vector r(0) (as weU, of course, as the 3-vector 
r (0)), namely 



f(0) = (a;(0) , 0,z(0)) , ^(0) = (0 , , 6) , 
entailing 

z{0) = a = f (0) • f (0) / f (0) = p/b , 

x{0) = c= r(0)Af(0) / f(0) =(a-a2)^^^ 



(A.21b) 

(A.21c) 
(A.21d) 



with a, P and b defined by (A.14d). This choice, using the representation (4.3) and (A. 19), 
also entails 



aFo{t) + bFi{t) cFo{t) -ibcF2{t) 
cFo{t)+ibcF2{t) -aFo{t)-bFi{t) 



(A.21e) 



{Beware: the quantities a and c introduced here, see (A.21a), (A21c), (A.21d) and 
(A.21e) have no relation to the "coupling constants" a and c in (A. 9b)). The standard 
matrix formula 



{d/dt) log {det [W (t)]} = Trace | {t) [W (t)]"^} 



(A.22) 



implies, see (A. 20b) and (A. 21a), that the determinant of the matrix W_{t) is time-inde- 
pendent. Hence from (A. 20c) we get 

det\W{t)] = l, (A.23) 

entailing that for W_ (t) and its inverse we can use the convenient parameterizations 



fit) -Ht) 

9 it) fit) 



fit) Ht) 

-9 it) fit) 



with 



f it) fit) + git) g it) 



(A.24) 



(A.25) 



(Beware: the scalar function f{t) introduced here should not be confused with the 
scalar/matrix function / introduced above, sec for instance (A. 2a)). 
Moreover, from (A.24) , (A.21) and (A.20) we get 



fit)=i[h{t)/2][ia + bt)f{t)+cg{t)] , 
g{t) = i[hit)/2] [- {a + bt)g{t) + cf{t)] 



(A.26a) 
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/ (t) = -i [h [t) /2] [(a + ht) f (t) + eg {t)\ , 

I it) = -i [h {t) /2] [- (a + bt) g it) + cf it)] , 
/(0) = /(0) = l, g{0)=g{0) = 0. 
Hence, if h (t) is real, as we hereafter assume, we see that 

fit) = rit), git)=g*it), 



(A.26b) 
(A.26c) 

(A.27) 



so that the remaining task is to solve (A. 26a), of course with the initial conditions (A. 26c). 
Once this task is achieved we get (from (A.20a), (A.21a), (A.21e) and (A.24) with (A.27)) 

Fo (t) = ^{-2 (a + bt) [f (t) g (i) + /* (t) g* (t)] 

+ c [f {t) + {t) - g' {t) - g*' it)] } , (A.28a) 
Fi it) = i{-aFo it) + {a + bt) [f (t) f* (t) - g (t) g* (t)] 

+ c[f{t)g*{t) + r{t)gm, (A.28b) 
F2 (t) = ^{2 (a + bt) [f it) g (t) - f* (t) g* (t)] 

- c [f it) - it) - / {t) + 5*' it)] } . (A.28c) 
To express these formulas in terms of / only, we moreover use the relations 

g{t)g*{t) = l-f{t)r{f) (A.29) 
(see (A. 25) and (A.27)), as well as 

g{t) = - [{a + bt) f{t)+i [2/h it)] f (t)] /c (A.30) 

(see the first of the (A.26a)), and its complex conjugate. We thus obtain the following 
expressions: 

1 



Foit) 



2c2/l2 (i) 



[h^ it) [c^ + {a + bt)'] [f it) + /*2 (t)] + 4 [p it) + /*2 (t)] } 



(A.31a) 



Fi{t) = l{-{a + bt) - aFo {t) +i[2/h (t)] [f (t) f* (t) - f* {t) f {t)] } , ( A.31b) 

(A.31c) 



These formulas provide explicit expression of the three functions Fj (t) in terms of the 
initial data (see (A.14a), (A.14b), (A.14d), (A.21c) and (A.21d)) and of the single function 
/ (t); they provide the solution to the initial value problem for (A. 13) via (A. 19). 

There finally remains to calculate / (t) . The most convenient route is by inserting 
the expression (A.30) in the second of the (A. 26a), obtaining thereby the following linear 
second- order ODE: 



f = Ih/h) f+ - {a + 2(3t + bh') {h/2y + ibh/2 



f 



(A.32) 
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with the initial conditions (see (A. 26c), (A. 30)): 

/(0) = 1, f{0)=iah{0)/2, 
or equivalently, via the position 

f{t) = [h{t)]'/\{t) , 
the foUowing hnear second-order ODE of (stationary) Schroedinger type: 



ip = ip 



ibh/2 - (a + 2pt + hH^) {h/2Y + - {h/h\ 



\ 2 



-h/h 



with the initial conditions: 

-1/2 



(^(0) = (/l(0))-^/^ <^(0) = (/i(0))~^/^ iah{{))-h{G) 12. 



(A.32b) 
(A.33) 

(A.34) 

(A.34b) 



The constants a ,/? and h are defined by (A.14d), and we have used (A.21c,d). The 
function hit) is of course defined by (A. 15) with (A. 14b), with /i (r^) defined by the 
original equation (A. 13) we set out to solve. 
For the simplest choice 



h (r^) = h{t) = k , 
(A. 32a) reads: 

/ = ^ [-k {a + 2pt + hH^) + 2ib] f , 

and the corresponding initial conditions (A. 32b) read 
/(0) = 1, f{0) = iak/2, 



(A.35a) 



(A.35b) 



(A.35c) 



It is amusing to note that (A. 35b) is, up to trivial changes (involving however the in- 
troduction of complex parameters), just the stationary Schroedinger equation for a one- 
dimensional harmonic oscillator, of course with the time variable playing the role of (one- 
dimensional) space variable. Its solution satisfying the initial conditions (A. 35c) can be 
easily expressed in terms of the parabolic cylinder function Dp (T)(see f.i. eq. 9.255 of [10]) 
and it reads: 



f{t) = ADp{T) + BDp{-T) , 



where: 



p = - (1 + ikc^/b) , 

T{t) = {ikhf^ {t-p/b'^) , 



A 



B 



D'J-t{0)) - kDj,{-T{0)) 



DpiTiO))D'p{-T{0)) - D,{-r{0))D>piTiO)) ' 

kDpi-TjO)) - D'pi-TjO)) 
Dp{T{0))D'^{-T{0)) - Dp{-rmD'^{T{0)) ' 



(A.36a) 

(A.36b) 
(A.36c) 

(A.36d) 
(A.36e) 
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k = ^ {ik/bf/'^ . (A.36f) 

Note that this case was studied in [3] where moreover it is mentioned that E.V. Fcrapontov 
already uncovered by geometric methods the relation with parabolic cylinder functions. 
The second simple choice we consider here is 

h [r') = A , (A.37a) 

entailing (see (A.14c), (A.15)) 

h{t)= k/[a + 2pt + bh'^] = k/x {t) ■ (A.37b) 
Then setting: 

t = -{2ibcz + p -ibc) /b'^ , (A.38a) 
f{t)=f{t{z)) = ]{z) , (A.38b) 

it is easily seen that (A.32a) reduces to the hypergeometric differential equation : 

z{l-z) f" (z) + (1 - 2z) f' (z) -k{k- 2ib) (26)^^ / (z) = , (A.38c) 

where the primes denote derivatives with respect to z. Thus the solution / (t) of (A. 32), 
with the initial conditions (A32b), can be easily obtained , via (A38a,b), in terms of the 
two independent solutions of (A38c) (see 9.151 and 9.153.2 of [10]). 
The next case worthy of consideration corresponds to the choice 

h{t)= k{a + 2pt + bh^) = kx~^''^ , (A.39) 

This case is significant from a physical point of view, since it is easily seen that, with this 
choice, (A. 13) becomes just the Newton equation of an electrical charge in the magnetic 
(radial) field of a magnetic monopole (or of a magnetic monopole in the electrical field of 
an electrical charge): 

f (t) = k[r (t)]"^ f{t) A r{t) . (A.40) 
Via the change of independent variable 

r= {a + 2pt + bHY' = [x{t)]'^' (A.41) 
(A. 32a) reads: 

(r^ - c^) /" (r) + r {3t^ - 2^) f' (r) - ^ (2z6r - k) f (r) = , (A.42) 

where of course / (r) = / (t) and primes denote derivatives with respect to r. An inves- 
tigation of this case is postponed to a subsequent paper. 

Let us finally report the solution of (A. 32), hence of (A12d), with 

h{t) = 2k/{l + 2iot) (A.43) 
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(a motivation for this choice will be clear in the next Appendix B, see (B.7d)). Let us set 



t = {2iz- bk) I (2 hkuj) . 
Then (A. 34a) becomes the Whittaker equation (see eq. 9.220.1 of [10]) 

1/1 \ \ 

z 



where the primes denote derivatives with respect to z and 

(p{z) = ip {t {z)) = (p{t), 

X = [26a;2 -ik {b^- 2/3 w)] / {Abu^) , 

ix''=e {4.f3u-b^ -Aau'') / {Alo^) . 



(A.44a) 



(A.44b) 



(A.44c) 
(A.44d) 
(A.44d) 



Thus, through (A.33), (A.43) and (A.44), the solution f (t) of (A.32) can be easily ex- 
pressed in terms of the two independent solutions of (A.44b), namely the Whittaker func- 
tions VFa.m (*)) and VF-a,m (--^ (0) (see eq. 9.220.4 of [10]). 



Appendix B. On the matrix evolution equation t/ = 2a^ + 



hU + c 



u, u 



In this appendix we discuss the linearizable matrix evolution equation (see (3.35)) 



U = 2aU + bU + c 



U, U 



(B.l) 



with a,b and c three arbitrary "coupling constants" (c could of course be rescaled away), 
and we identify two cases in which, at least for 2 x 2matrices, this equation can be explicitly 
solved in terms of known special functions. We set 

Uit) = 'Pit) V{t) , T = T{t) , (B.2) 

and thereby obtain 

y" = V' [2af - 2f (if/^) - r] /f^ 

+ V[b + '2ai<p/^)-{^/<p)]/f^ + c{^/f)[v:,V] . (B.3) 

Here and throughout this appendix primes denote derivatives with respect to r. 
A choice naturally suggested by this equation is 

f{t) = ip{t) , (B.4a) 

which entails that F (r) satisfies the matrix ODE 

V' = ip'^[2aip-3ip]V! + ip'^ [2aip + bip-ip]V + c[V',V] , (B.4b) 
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which we report here for future memory. However, we prefer now to focus on two special 
cases of (B.l), which allow a completely explicit solution (at least in the case of (2 x 2)- 
matrices). Hence we set 

if (t) = exp (fit) , (B.5a) 
f{t) = exp(i/t) , (B.5b) 

so that (B.3) becomes 

V" = p exp {-ut) V + gexp (-2 i/t) F + cexp [(/x - v) t\ [F', Y\ , (B.5c) 

with 

p = 2{a- fi) - , (B.5d) 

g = 6 + 2o/x-jU^ . (B.5e) 

The first case we consider is characterized by the condition 

b = - (8/9) (B.6a) 

to which we associate the choices ^ = v = 2a/3, which entail, via (B.5d), p = and, via 
(B.5e) and (B.6a), q = Q as well. Hence in this case the position 

£/(t) = exp(2at/3) Z(r) , (B.6b) 

T {t) = [exp (2 a t/3) - 1] / (2 a/3) , (B.6c) 

yields 

Z" = c [Y!, V] , (B.6d) 

which, up to a notational change and a trivial rescaling of the dependent matrix variable, 
is just the matrix evolution equation discussed in Appendix A [solvable in the 2 x 2case, 
see (A.12d), (A. 35a), (4.3), (A. 19), (A. 28) and (A. 36) ). Hence a convenient prescription 
to solve (B.l) with (B.6a) is to solve instead (B.6d), and then perform the change of 
dependent and independent variables (B.6b) with (B.6c). Clearly, this implies that all 
solutions of (B.l) are completely periodic with period T = 27r/u; if 

a = (3/2) iu, h = 2j^ , (B.6e) 

with u an arbitrary (nonvanishing) real constant. Note however that in this case the 
matrix evolution equation (B.l) is complex. 

The second case we consider is characterized by the restriction 

6 = 0. (B.7a) 

In this case we set = 0, = 2a, which again entails, via (B.5d), p = and, via (B.5e) 
and (B.6a), = as well. Hence in this case the position 

U{t) = V{T) , (B.7b) 
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T(t) = [exp(2at) -l]/(2a) , (B.7c) 

yields 

y" = c (1 + 2 ar)-^ [Y!, V] , (B.7d) 

which, up to a notational change, a trivial rescaling of the dependent matrix variable, 
and a redefinition of the independent "time" variable, is the solvable matrix evolution 
equation discussed at the end of Appendix A (see (A.12d), (A.43), (4.3), (A.19), (A.28) 
and (A. 44)). Hence a convenient prescription to solve (B.l) with (B.7a) is (at least in 
the case of (2 x 2)-matrices), to solve instead (B.7d), and then perform the change of 
dependent and independent variables (B.7b) with (B.7c). Clearly, this again implies that 
all solutions of (B.l) are completely periodic with period T = 2-11/10 if 

a = iio/2, 6 = 0. (B.7e) 

Of course in this case as well the matrix evolution equation (B.l) is complex. 

Appendix C. Proof of a formula 

In this appendix we prove a formula relevant in connection with the first of the two 



techniques introduced at the end of Section 3.5 to get, from an evolution equation involving 



one, or a few, functions, (coupled) evolution equations involving many more functions 
(which might be matrix- valued) . 

Let / (z) be an analytic function of 2: , so that, at least for small enough \z\ , 

00 

m=0 

is well defined. 

Let r]j, J = 1, J be J elements of an Abelian algebra such that 

rjjm = VkVj = Vj+k, VJ = 1, r]j±J = Vj ■ (C.2) 

The last of these equations is consistent with the assumption made throughout this 
appendix, that all indices are always defined mod (J) . 
Then clearly the formula 

J 

= ^Vj fjizi,...,zj) (C.3) 
i=i 

defines uniquely the J functions fj(zi,...,zj) . The purpose of this appendix is to derive 
an explicit formula to compute these functions. But before doing so, let us emphasize that 
(C.3) and (C.2) entail (with self-evident notation) 

\fcl=l / \k2 = l J j = i k=l 

(C.4a) 
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\ki=l J \k2=l J \k3=l 

J J J 

= X ''i X X fj-k-k' {zi, zj) Qk {zi, zj) hk' {zi, zj) , (C.4b) 
j=l k=l k'=l 

and so on. 

To compute the functions fj{zi,...,zj) , see (C.3), we note that the algebra (C.2) 
admits the following J realizations: 

=e^p{27rijK/J), j = l,...,J, (C.5) 
with K = 1, J. These realizations are all different iff J is prime. Hence (C.3) entails 

\A;=1 / j=l 

namely 

/ (Y.Zke^p{27rikK/J)\ = ^ /j' -> ^j) exp (27^^i if/J) . (C.6b) 
\ifc=i / j=i 

Multiplication by exp {—2 -rr i k K / J) yields 



exp{-2TrikK/J) / (^E Zk'e^p{2Trik'K/j)\ 

\k'=l ) 



J 

X fj -,^j)exp [27rz (i - k) K/J] . (C.7) 



We sum now over K from 1 to J , and use the formula 

J 



J2(^^v{'27rijK/J) = JSj,j, j = l,...,J, (C.8) 



K=l 

getting thereby 



fj{zi,...,zj) = J-^ X exp{-^^ijK/J)f(^Zke^p{27rikK/J) \ , (C.9) 

K=l \k=l J 

which is the formula we intended to get. 

The generalization to the case when / is matrix- valued is obvious, as well as to the 
case when it depends on more than one argument, say / = / (x; y); then 

(J J \ J 

^ miXkil^VkiVki] = ^Vj fj{xi,-,xj;yi,...,yj) (C.IO) 
ki=l k2 I j=i 
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defines the J functions fj (xi, ...,xj;yi, ...,yj), and 



fj{xi,...,xj;yi,...,yj) = J ^ ^ exp (-2 7ri j K/ J) 



K=l 



f ( ^ Xk^exp{2'KikiK/ J); ^ y^jexp (27r i A;2 -fC/J) 

ifci=l A;2=l 



(C.ll) 



/O 1 
^-(,10 

- i {r ■ a)] 



I i 

^y=\-^ 



( 1 

-1 



Appendix D. Algebra of a-matrices, and a useful 3-vector 
identity 

In this appendix we display some standard formulas for the a-matrices, and a useful 3- 
vector identity. 

(D.l) 

(D.2) 
(D.3) 

(D.4) 

(D.5) 

(D.6) 
(D.7) 

(D.8) 

(D.9) 
(D.IO) 

(D.lla) 



[p + i {r ■ a}] ^ = [p — i {f ■ q_ 

^(r-<i) ■ a) , (r-<2) ■a)]=2i (f^^^ A 7<^^) ■ a , 
(r-<^) • a) (r-<2) • a) (r^'^ ■ a) = i (f^'^ A f^'^) ■ f^'^ 

' ' ' '"^ ' ' f(3) j^^i) .fK2)y 



= 2 [r<^) fr<2) .r<3)^ 



(P + 



\ r.-<3) . ^ + (^^-<3) . ^ (^2) . j^^-<i) . ^ 

I _ ^2) j^^l) . ^3)^ ^ ^3) . ^2)^ 

. f(2) j _ f(2) 



r = —if 



- i r 



^1) . ^3)' 

2r-<i) 

.)i/2_(^_i^)i/2]/(2r) , 

-\- i (f- ?t\ sin ( r\ . 



if-q)^^'^ = p + if-ff_, p= {p + ir)^^^ + {p — ir)^^^ /2 



'jj •£ , 

(r-<i)T-<^))" 



exp {if - a) = cos (r) + i (r • a) r ^ sin (r) 
/ (r • a = {/ (r) + / (-r) + [/ (r) - / (-r)] f- a/r} /2 

The 3-vector formula 

ar + r A 6 + • (/| c = f 

entails 

f=[a^ + b^y^ af + bAf + jbAc + ajc + 
c?f- d + a(^f Actj -1)+ (^f-b^ • 5 



7 = 



(^f-b + jc-tj b , (D.llb) 
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(a^ + 6^) + a^c- d + a (c/\ dj ■ b+ (c- 6) (d^- 6) 



-1 



(D.llc) 



The subcases 6 = , c = (or d = 0), a = 0, while easily obtainable from the above 
formulas, deserve separate display. Case 6 = 0: 



r = a 



f-[f-d)(a + c-d) 



Case c = (or d = 0) : 

r=[a^ + b^)~^ af + bAf + a-^ (^f-tj 

Case a = (note that it requires 6 ■ c ^ as well as 6 • d ^ ) : 

-1 ^ -1 



(D.lld) 



(D.lle) 



'{bAf- (b-c 



[b-pj bAc + (b-c 



{b.dj 



62 (6 • /) - (b-c) (bAf) ■ d+ (b-f) (bA c) • d] 6} 



(D.llf) 



Appendix E. Some solvable and/or linearizable matrix evo- 
lution equations 

In this appendix we report several solvable and/or linearizable matrix evolution equations. 
We also outline their derivation, although none of these results is really quite new (see, 
for instance, [11]). Most of these equations are of "nearest-neighbour" type; as always in 
this paper we ignore here the question of the boundary conditions to be assigned at the 
extremal values of n (say, for n = and n = N + 1). 

Let us take as starting point the following solvable linear matrix equation with constant 
(time- independent, matrix) coefficients: 



W„ (t) = lEn (t) + B„ W„^, {t) + W^_, (t) , 
and let us set 



—n -LL-n+1 -!-!-n 



(E.l) 



(E.2) 



(here and below wc often omit, for notational simplicity, the explicit indication of the time 
dependence) . Then the matrix (t) evolves according to the nonlinear equation 



Vn {t) = Zn {t) - Zn (*) An + (0 Z„ (t) 

- Vn it) it) + Cn+l - Zn (*) [Zn-1 (*)] 



(E.3) 



The special case of this equation with A^^ = A, B^^ = 5 and C„ = is the first 
nontrivial evolution equation of the so-called discrete Burger's hierarchy [11]. 

There are now various ways to derive, from this first-order solvable matrix evolution 
equation, solvable or linearizable second-order matrix evolution equations. We describe 
two of them. 
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A first trick is to separate the odd/even labeled matrices, by setting, say, 1^2m ~ 
Z2m+i = U.m^ A2m = Am^ ^2m+i = lm> and SO On. This yields 



Hn = AnUn- Hn An + B^Hn Hn ' Hn Hn + Qr, - C„ 



Hn-l 



(E.4a) 



= in+l U.n -RnAn+ R^+l Hn+l U.n ' U.n U.n + - Rn Q-^ \^n\ 



-1 



(E.4b) 



We then time-differentiate the first of these two equations, use the second to eliminate 

l}_ , and use the first (undifferentiated) to eliminate IJ_ after having set, for simplicity's 
sake, C„ = (the diligent reader is welcome to work out the more general result that is 
obtained forsaking this simplification). We thus get: 



An+l + lLn+1 



B. 



-1 r 



ILn ~ AnU.n~^ ILn An + ILn ^n H-n ~ Q-r, 



Un -AnUn+ UnAn + UnRnUn ' 



An + B^Hn 



(E.5) 



Let us now consider the simpler case in which the constant matrices are replaced by 
scalars, namely we set = an 1, An = ^nl > and so on. We thus get 



Un = [(^n - an+l] Cn + [a„+l - On] [ttn - a„] + C„ [bnUn " bn+llLn+l] 
+ [Un - dn] bn+1 ILn+l U-n - '^Ktln Un 
+ [a„ + a„ — 2 On + hn+i Un+l-bnUn] £n + I^n • 



(E.6) 



Let us also record the special case of this solvable matrix evolution equation which 
corresponds to the position (consistent with (E.6) ) C/„ = IL, Cn = —c/a, bn = b, an = 
d + an, dn = a — a + an. This yields 



U = c + aaU+{a-a)U + baU^-2bUU, 



(E.7) 



a solvable evolution equation which contains the 4 arbitrary constants c, a, a, b. 

A second method to obtain a second-order equation from the solvable first-order equa- 
tion (E.3) is by setting 



Vn it) = M„ it) [M„ (t)]-' , M„ it) = Zn it) M„ it) . 
One obtains thereby the following evolution equation for the matrix (i) : 

Kn = Kn Kn + An+l Kn ' Kn [Mn]~' AnK^ 

+ B^+lM^+l [M„+l]"' Kn-Kn [MJ"' S„M„ 



(E.8) 



M 



■n—l 



(E.9) 



This should be categorized as a linearizable system of matrix evolution equations, since 
to solve it one must solve, in addition to the solvable equation (E.l) (to get it) and 
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then, via (E.2), (t)), a linear nonautonomous matrix evolution equation (to get M „ (t); 
see the second of the (E.8)). 

Let us call attention to the presence, in this equation, of a time-differentiated matrix 
in the denominator (see the last term in the right-hand-side). 

As an obvious consequence of the way this equation, (E.9), has been derived, see (E.8), it 
is invariant under the transformation M „ — > M^D^ with D_n arbitrary constant matrices, 

= 0. 

When all the constant matrices are replaced by scalars, namely if we set in (E.9) 
= a^l, = bnX: Q.n = we obtain the following linearizable matrix evolution 
equation: 



M„ = K+i - an - bn] Kn + {Ku [KJ ' + M„+i [M„+i] M„ 



+ Cn+1 



M. 



.n-l 



-1 



(E.IO) 



Let us end by noting that the (compatible) position M„(t) = M_{t), On = n (2a + c— 1), 
bn = c — 1, Cn = nb yields for the matrix M_{t) precisely the solvable evolution equation 
(3.1). 



Appendix F. A linearizable matrix system, its solvable ver- 
sions and a class of solvable Hamiltonian many-body prob- 
lems 



In this appendix we show that the system of matrix evolution equations 

N 



{(^nm l£.ni 



N 



UnJ{Uj,v^;t) 

VnJ{Uj,V,;t) 



(F.la) 
(F.lb) 



is linearizable. Here the 4 A'^^ quantities anmi ^nm: Cnm: dnm are arbitrary (they could also 
be time-dependent functions), and F is an arbitrary function of the 2 matrices C/^, V-m 
and of the time t; note however that the same / enters in (Fla) and (F.lb), that this quan- 
tity is independent of the index n, and that it is a scalar/matrix function of its arguments, 
namely it satisfies the property W^/ ([7^,^^; i) W_-^ = l{WUjW-\ WVjWr^^t). 

To prove that (F.l) is linearizable we proceed as in Section 3.2, namely we set (see 
(3.11) and (3.13); beware of the notational changes!) 



w{t)Un{t)mt)r\ Unit) = mt)r' Unit) w{t) , 

Vnit)=W{t)Vnit)[W{t)r\ Vnit) = \Wit)r' ILnit)W{t) , 

W{t)=W (t) / [[/, (t) , F, (t) ; t] , W{t)=l [n, (t) , vj (t) ; t] W {t) . 
Time-differentiation of the first of the (F.2) yields (using the first of the (F.3)) 



iLn = W {Un - [Un, / iHj, t)] } WT^ , 



(F.2a) 
(F.2b) 
(F.3) 

(F.4a) 
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w{v^-[v„J{Uj,Vj;t)]]w- 



(F.4b) 



and from these equations and (F.l) , (F2) we see that the matrices u„ (t) , u„ (t) satisfy 
the linear evolution equations 



TV 



N 



(F.5) 



m=l 



m=l 



The linearizabiUty of (F.l) is thereby proven, since its solution can be achieved by solving 
firstly (F.5), then the second of the (F.3) (a linear equation for the matrix W_{t)), and 
then recovering C^„(t), y_^{t) from the second of the (F.2a,b). 

It is moreover easily seen that in some cases, see below, the system (F.l) is not only 
linearizable, it is in fact solvable. An obvious (and rather trivial) case is if the quantities 
anm, Km, Cnm, dnm all vanish, a„m = hnm = Cnm = dnm = 0, and / does not depend 
explicitly on the time t, f {V_j, V_j;t) = f (U_j, F^). In such a case (F.5) entails that the 
matrices Un (0 > Hn (0 f^'^^ time-independent, u^(t) = li„(0), u„(t) = ?Z„(0), and the 

second of the (F.3) becomes explicitly solvable. Hence one concludes that the equations 



U„ = 



v„ = 



are explicitly solvable: 

{t) = exp [-tl{Uj (0) , (0))] (0) exp [tf {U^ (0) , (0)) 

(t) = exp \-tl ([/, (0) , y, (0))j Zn (0) exp / ([/, (0) , F, (0))' 



(F.6a) 

(F.6b) 
(F.6c) 



Another, perhaps less trivial, case in which the equations (F.l) are in fact also solvable 
obtains if 



N 



(F.7a) 



and the quantities a„^, bnm, Cnm, dnm, are all time-independent and satisfy the con- 
straints 



(F.7b) 



Indeed it is easily seen that these conditions arc sufficient to guarantee, via (F.l), that 
/ is time-independent, / = 0, so that there holds again the explicit solution (F.6b,c), of 
course with / given by (F.7a). An interesting case (see below) is that with A„ = A, so 
that dnm = —o-mn and (F.l) read 



N 



Un=J2 i^nmUm + bnmVm + ^ lUn, Um^ Z^]]) , 
m=l 
N 

Vn=^ {CnmUm " Omn Zm + ^ [Zn> \^m^ Y-mW) ' 



(F.7c) 



m=l 
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with 

(F.7d) 



-'ran ; ^nm 



Linearizable and solvable evolution equations in ordinary (3-dimensional) space can be 
obtained from the above matrix equations, by using (compatible) parameterizations in 
terms of 3-vectors, see Section We only display the equations that are obtained from 
the solvable equation (F.7c) with (F.7d) via the parameterization (4.3), namely by setting 
= iqn ■ V_^ = ipn • (7- They read: 

N 



Pm + 4:X [g„ A [qm A Pm]]) , 



m=l 
N 

P„ = ^ (Cnm Qm " Omn Pm + 4 A [pn A [qm A Pm]]) ■ (F.8a) 
m=l 

Remarkably, these solvable equations of motion, with (F.7d), are precisely the equations 
of motion entailed by the Hamiltonian 

N 

^ iQjiPj) — ^ ^ {o-nm Qn ' Pm ~l~ bnmPn ' Pm/^ Cnm Qn ' Qm/'^ 
n,m=l 

+ 4 A (gn A Pn) ■ {qm A Pm)} ■ (F.8b) 

Note that the {2N^ + N + l) constants that appear in this Hamiltonian, and in the 
equations of motion (F.8a), are arbitrary. But the solvability of this Hamiltonian system 
has a rather trivial origin; indeed, the evolution equations (F.Sa) are hardly nonlinear, 
since they entail, as it can be easily verified by direct computation (and it is of course 

N 

implied by our treatment), that the quantity ^ {qm f\Pm) is a constant of the motion. 

m=l 

Nevertheless, even for = 1, when the Hamiltonian reads simply 

H{q,p) = aq-p + bp'^/2-cq'^/2 + 2X{qAp)'^ , (F.9a) 

and the corresponding equations read simply 

q = aq + bp + 4X [qA[qAp\], p=cq-ap + 4X[pA[qAp\] , (F.9b) 

this is a rather interesting system, whose detailed analysis, however, exceeds the scope of 
the present paper. 

We end this appendix by reporting the many-body generalization of the last two equa- 
tions, (F.9b), that are obtained by using the first of the two multiplication techniques 
described at the end of Section (see also Appendix C). Hence we set 

N N 

P=^r]nPn, q = '^Vnqn, (F-IO) 

n=l 71=1 
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with analogous formulas for the "coupling constants" a, b, c, A (see (F.9)), as well as 

N 

H{p,q} = '^VnHn{pi,...,PN;qi,-,qN) , (F.lla) 



n=l 



with (see (F.9a)) 



N 

Hn{Pl,-,PN;qi,-;qN) = ^ [an-m-UiiPni ■ qn2j 

7ii,n2=l 

~l~ bn—ni—n2 {Pni ' Pn2)/'^ ~l~ Cn—n\—n2 {qni ' 9n2)/^] 
N 

+ 2 J2 ^" 

— ni— n2— na— 714 (Pni f\qn2)- (Pns ^ qrn) ■ (F.llb) 

ni,n2,n3,n4=l 

Here and below all the indices n, ni, are defined mod (N). 

It is then easily seen that the Hamiltonians are in involution, and that the 
Hamiltonian H]\r yields the solvable N-body equations of motion entailed via (F.IO) from 
(F.9b), namely 



N N 



q-n — ^ ^ (o-n— ni Q'rti ~l~ ^n— ni Pni) ~\~ 4 ^ ^ [An— m— n2— ns (Prii ^ '?«2) ^ ' 
ni=l ni,n2,n3=l 

(F.12a) 

AT TV 

Pn — ^ ^ (Cn— rii Q'ni ^n—n\ Pn\) ~\~ 4 ^ ^ [An— ni— n2— na {.Pni ^ qn2^ ^ Pn-s\ • 
ni=l ni,n2,n3=l 

(F.12b) 



But of course, as emphasized in Section |3.5| , these coupled equations of motion can be 
decoupled by linear transformations of type (3.92). 

Appendix G. Example: first-order and second-order lineariz- 
able matrix ODEs, and a linearizable one-body problem 

In this appendix we illustrate, via a simple example, the possibility, in the context of the 



technique of Section 3.2, to restrict firstly attention to first-order matrix ODEs and to 
obtain subsequently second-order ODEs by appropriate additional steps. We only con- 
sider an illustrative, very simple, example; this allows, at very little cost in terms of 
repetitiveness, a completely self-contained presentation, which can be understood without 



any previous knowledge of the treatment of Section 3.2; but for the diligent reader inter- 



ested in the connection we note that the case considered in this appendix corresponds. 



up to a trivial notational change {ui u, U2 ~^ IL)^ to the treatment of Section 3.2 with 
N = 2,Ui =ui-{aui + (5u2) , Uji = M2-(7«i +(^^2) > Z = /-(a^i + 6«2 + c [ui, U2]) 
and g = h = v = y = 0. At the end we display the (highly nonlinear) linearizable one-body 
problem which corresponds to the second-order linearizable matrix ODE via the parame- 
terization (4.3). 
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Let us set 



u = WUW-^, U = W-^uW, 
v = WVW-^, V = W-^vW, 

W = W {aU + bV + c\U, V]} = {au + bv + c \u, v\} W , 
u = au + Pv, v = ^u + Sv, 

with a, b, c, a, 7, 6 arbitrary constants. 

This entails for the matrices C/ and V_ the first-order nonlinear ODEs 

U = aU + pV + b\U,Y] + c\U, \U,V]] , 
V = jU + SV + a[V,U]+c[V,[V,U]] . 



(G.la) 
(G.lb) 
(G.2) 
(G.3) 



(G.4a) 
(G.4b) 



These nonlinear matrix ODEs are of course linearizable, since their solution can be 
achieved via the following steps: (i) set (for simplicity) W (0) = 1 as initial condition to 
complement (G.2); (ii) note that this entails u(0) = U{0), v{0) = V{0) (see (G.la,b)); 
(iii) evaluate u(t) and v{t) from the (explicitly solvable) evolution equation (G.3), taking 
into account the appropriate initial conditions, see (ii); (iv) evaluate W_{t) by solving 
the second of the (G.2), with initial condition 15^(0) = 1, see (i) (note that this is a 
linear nonautonomous matrix ODE, entailing the solutions of Manalogous systems of M 
linear first-order coupled nonautonomous ODEs - assuming we are dealing with (M x M)- 
matrices); (v) finally evaluate U_{t) and Y_{t) from the second of the (G.la,b). 

Let us now derive, from the 2 first-order ODEs (G.4) satisfied by the 2 matrices C/(t) 
and Y-{t)i a single second-order ODE for one of these two matrices, say for lJ_{t). This 
is easily obtained by time-differentiating (G.4a), thereby obtaining (using (G.4b)) the 
second- order linearizable ODE 

U = aU-\-PjU-\-/3SV+ (bS-aP) \U,V]-\-b \u, v] + {c6 - aP) \U, \U, V]] 



+ c 



u,v 



+ 



U, [U, Z]] - P m, V],V]}-c{a [U, [U, [U, V]]] 



+ b \u, [\u, Z] , Z]]} - \u, [u, [\u, Y] , m 



(G.5a) 



where the matrix V_ should be expressed in terms of IJ_ and U_ by solving for V_ the (non- 
differential) linear matrix equation (G.4a), namely 



PV-b[V,U] + c[[V,U],U] = U-aU. 



(G.5b) 



We end this appendix by displaying the linearizable one-body problem that corresponds, 
via the parameterization (4.3) (and by using (D.ll)), to the linearizable matrix ODE (G.5): 



f=ar-\-P'yr-\-PSr-\-2{Pa-b5) rAr-\-2brAr 

-\- 4{cd - afi) fA (^fA -|- 4c | fA (^f A -|- fA 

-|-4c{/3 (f A f) A f -I- 2ac fA^fArj Af 



-1-26 f A 



f A r A f 



16c^ f A 



r A f A f A f 



^f A I 



(G.6a) 
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where 



r = 



(/3-4cr^)^ + 46V ^ {2b f Ar+ {/3 - 4cr^) 'r 



2\ y 



+ 



46^ (/3-4c5 



-1 V 



r + (4c X — a) r ) 



and 



X 



+ (/? - 4c r^) (4c X - a)] f) 

(/3-4cr2)^ + 86cr2 (/?- 4c r^)^ + 46^/3 
(/3 - 4cr^)^ + 45^r^ (^a - f • . 



(G.6b) 



n -1 



(G.6c) 



Appendix H. Solvable matrix ODEs, and a solvable many- 
body problem 

In this appendix we manufacture a set of solvable matrix evolution equations, and we 
obtain the corresponding solvable equations of motion for 3-vectors by applying to the 
matrices the parameterization (4.3). 
We start from the equations 



N 



m=l 
LLn 



Here, for simplicity, we assume the N(3N + 1) quantities a 



nrm "nmj '-■nm 



(H.la) 

(H.lb) 

, an to be time- 
independent "coupling constants," although the more general case in which they are given 
functions of the time t could be easily treated as well. 

Firstly, let us obtain the second- order evolution equations satisfied by the N matrices 
ILn = U-n (^)' these equations are solvable, because we show next that the equations (H.l) 
can be solved by performing only algebraic operations. 
Time-differentiation of (H.la) yields, using (H.lb), 



m=l 



^nm Ujm. "I" 



U V 



~^ Cnm (^r, 



H-mi — m 



where we have introduced the convenient notation 

^nm — Oinm 0!m ~^ bnm ■ 



(H.2a) 



(H.2b) 



In (H.2a), the matrices j/„ are supposed to be expressed in terms of IJ_^ and L[„j by solving 
the (algebraic) equations (H.la), not the (differential) equations (H.lb). This can always 
be done algebraically (up to obvious restrictions, see below), but hereafter we restrict 
attention to the simpler case characterized by the restriction 



— dnm. OjTj 



Cnm — dnm Cm , 



(H.3a) 
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which expresses the 2N'^ constants a„^, Cnm in terms of the (iV + 2) (arbitrary) con- 
stants dnrm o-n, (hi- We moreover assume that the (AT x A?^)-matrix D, with matrix elements 
dnm, is invertible, and we term dnm the matrix elements of the inverse matrix D = D: 

dnm = (d) , dnm = (r') = {D)^^ . (H.3b) 

Then clearly from (H.la) we get 



N 



m=l 



N 



m'=l 



(H.4) 



a linear (nondiffcrential!) equation for the matrices y„ which can be conveniently solved 
in explicit form [6] (see below). 

Let us now show how to solve (H.2), or rather, equivalently, (H.l). From (H.lb) we get 

Vn{t) = anUnit) + Cn, (H.5a) 
with the constant matrices C„ given in terms of the initial data: 

C, = Z„(0)-anf/n(0) . (H.5b) 

Insertion of (H.5a) in (H.la) yields the following set of linear ODEs with constant 
coefficients for the matrices U_n{t), which can of course be explicitly solved by purely 
algebraic operations: 

N 

Um + CJ + bnm Hm + Cnm , Cj} . (H.6) 

m=l 

Next, let us exhibit the solvable equations of motions which are obtained by applying 
the parameterization (4.3) to the solvable matrix evolution equations (H.2a) with (H.4) 
(we leave as an exercise for the diligent reader the derivation of the more general equations 
that obtain via the parameterization (4.1) rather than (4.3)). First of all we note that 
(H.4) can now be rewritten as follows: 

On Vn +7nVnAfn = Wn , (H.7a) 

where we set, for notational convenience, 

7n = 2c„, (H.7b) 
as well as 



N 

Wn = 

m=l 



N 
m'=l 



(H.7c) 



and of course C/„ = ifn ■ ff, F„ = ivn ■ Hence, using (D.ll), 

= [7n (r„ • Wn) + ItJn + On In Tn A Wn] / [an {an + 7n '^n)] • (H.7d) 
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On the other hand, via the parameterization (4.3), the solvable matrix evolution equations 
(H.2a) read 



N 



m=l 



(H.8) 



where we used (H.3a) and (H.7b). Here the N three- vectors t/nshould be expressed via 
(H.7d,c). 

Let us write in more explicit form these equations of motion in the special case in which 
the matrix is diagonal, indeed when we set (without any additional loss of generality) 



D = D = h d, 



nm — "nm 



(H.9) 



We thus get the following solvable many-body problem, which only features two-body 
(velocity-dependent) forces, and contains the N (N + 3) coupling constants a„, a„, 7„, 
bnm, which are arbitrary except for the restrictions a„ 7^ (see the following equations of 
motion) : 



Tn = anfn + an 7n r„ A "t" {-Jn Tn A 7„ [fn ■ fnj Vn + a„ 7n r„ A r„ 
N 

^ ^ {bnm {(^m fm Tn ^ \ln (Tn ' ^m) "I" ^m "I" Tn ^ ^m] }}}/ 



m=l 



[an K + Tn rn)] ■ 



(H.IO) 
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